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ABSTRACT 



Many of the advances in servo-system design fostered by the general 
availability of digital computers are not utilized by practicing engineers 
because they discard rather than complement existing methods. An exception 
is the parameter plane method which retains the principle of design by such 
familiar specifications as damping and settling time, but where classical 
methods display the effect of only one variable on these specifications the 
parameter plane method can display the simultaneous effect of two. In this 
paper the restrictions on problems suitable for parameter plane analysis 
are eased and the number of specifications which can be considered on the 
parameter plane is enlarged to include bandwidth and steady state error. 
Parameter plane methods are also adapted to dominant mode design in a 
manner which not only allows selection of the dominant mode but simultan- 
eously guarantees its dominance. These Extensions of parameter plane theory 
are related to both continuous and sampled data systems. 
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CHAPTER I 



In recent years great emphasis has been placed on the "controls field" 
in electrical engineering. With this emphasis has come a proliferation of 
techniques for analysis and design of servomechanisms, many of which have 
not yet been fully exploited. Among such techniques are the parameter plane 
methods which can be extended in both applicability and in specifications 
considered. 

Before analyzing the parameter plane as a design tool it may be wise to 
establish its place among tools presently available to the engineer. A 
broadbrush review of available methods should be sufficient for this pur- 
pose. Present day servomechanism design subdivides conveniently into two 
disciplines, direct time domain analysis and frequency or transform analysis. 
Of the two the former appears to be the more powerful albeit the more com- 
plex. Design by this procedure ordinarily includes the steps: 

a) determination of signal characteristics 

b) selection of a criterion by which system performance can be 
evaluated (i.e. integral squared error, etc.) 

When the results of these steps are clear cut and definite time domain 
analysis will probably be the best tool available to the design engineer. 

This holds particularly for sampled data systems where incorporation of a 
digital computer can make realizability of theoretically obtained results a 
minor problem. Despite its power time domain analysis is by no means a 
panacea. By its very nature the time domain approach demands that system 
structure, its inputs and its objectives be well defined. These are by no 
means trivial requirements since a system designed on the basis of integral 
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squared error, for example, need not have acceptable acceleration character- 
istics, overshoot, frequency response, etc. Furthermore, analysis in the 
time domain generally provides little insight into possible performance 
criteria trade-offs or structural changes which might enhance the overall 
system. 

The difficulties of obtaining performance characteristics from the 
integrodifferential equations which describe a servo-system have to a large 
extent been overcome by the use of transforms which reduce the operations 
of differentiation and integration to operations of multiplication and 
division. When transformed, the pole zero pattern or magnitude phase plot 
of the integrodifferential equation furnishes an indication of system per- 
formance by inspection. In addition these plots often indicate system 
modifications or compensation which will improve system performance. The 
primary disadvantage of transform analysis is that pole-zero patterns, 
bandwidth and other transform criteria abstract the time domain performance 
of the system and it is in the time domain that the system must ultimately 
perform satisfactorily. 

A great number of servo design procedures fall within the category of 
transform techniques. In a broad sense these procedures can be subdivided 
into a class dominated by a synthesis concept and a class dominated by a 
compensation concept although in practice neither can be isolated from the 
other. The synthesis approach to servo design is accomplished in three 
steps. 

1) determination of the closed loop transfer function from 
specifications 

2) determination of the open loop transfer function from the 
closed loop transfer function obtained in step (1) 



2 









-•I 







3) synthesis of the open loop transfer function with due regard 
to components which are not at the engineer's disposal. 

While the transformed equations are sufficient to determine the system's 
time domain behavior the complex relationship between the two thwarts ready 
selection of the "best" transfer function. Hence it is likely that the 
selected transfer function which meets specifications is more difficult to 
synthesize than some other transfer function which also meets specifications. 
Furthermore, there is no systematic or iterative procedure which will im- 
prove this initially selected transfer function. To this liability must 
be added the problems that fixed system components impose on the choice of 
possible closed loop transfer functions. To its credit the synthesis con- 
cept follows a logical sequence from specification to final realization. 

This is in contrast to the compensation approach which attempts to "patch 
up" an unsatisfactory initial system in the hope it can eventually be made 
to meet specifications. The second liability has also been partially elimin- 
ated at the expense of simplicity. With the parameters at his disposal the 
engineer is able to produce the system which most closely approximates the 
one Initially chosen. [1] However this closest point of approach is not 
guaranteed to meet specifications nor does it preclude that a satisfactory 
system could be obtained with available parameters were a different choice 
for the "ideal" closed loop transfer function made. 

Utilization of transform techniques to design compensation for an 
existing skeleton servo-system is an approach more widely employed than syn- 
thesis. Two steps are basic to design by compensation. 

1) From specifications a "best choice" for the fixed elements of 
the servo-system is made 
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2) Then compensation for these elements not under the engineer *s 
control is designed to insure system compliance with specifi- 
cations* 

These steps contrast vividly with those employed in synthesis, largely elim- 
inate the drawbacks of that method and offer very definite advantages. The 
realization problem is eliminated by the expedient of employing only realiz- 
able compensators to modify the fixed skeleton system. Of almost equal im- 
portance are the well established methods available for design of compensa- 
tion* They are for the most part simple, powerful and provide insight in 
addition to numerical solutions; a feature which is almost unique to these 
methods* In addition systems designed from a compensation concept tend to 
be comparatively inexpensive, a fortuitous result which usually occurs when 
dealing with the physical compensation rather than the theoretical system. 
However design by compensation also has its disadvantages* Most serious is 
ignorance of whether the possibility for compensation exists; the fixed 
system may be incapable of being compensated to meet specifications. In 
this respect the s)mthesis concept is clearly superior for whenever an accept- 
able theoretical system cannot be found then it serves no purpose to attempt 
design of an impossible physical system* 

The preceding thumbnail sketches suffice to acquaint the reader with 
approaches available to the servo-system designer* However, deficiencies 
inherent to an approach are subject to the **state of the art** and to the 
state of development of the methods employed* The compensation approach 
illustrates the case in point. In an abstract formulation its merits are 
not outstanding yet the simplicity of the techniques and the well documented 
correlations between available variables and design criteria serve to make 
this the most prevalent of approaches. Parameter plane techniques also 
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advantageously utilize these well established correlations between f re quency 
response, bandwidth, pole -zero patterns and time domain performance specif ica 
tions such as rise time, maximum overshoot and settling time. Furthermore 
the essential steady state nature of frequency response analysis and the 
one parameter limitation of root locus diagrams are largely overcome by 
parameter plane methods. At the expense of simplicity the effect of simultan 
eous variation of two parameters on root locations can be displayed. Band- 
width specifications and steady state error requirements can also be indicat- 
ed on this plot. Furthermore such troublesome questions as second order 
dominance or whether dominance is possible by adjustment of two parameters 
can be answered by inspection. These advantages auger more widespread applic 
tion of parameter plane methods in systems of up to moderate complexity, 
especially as the digital computer becomes more readily available to compute 
the required curves. 

The intent of this chapter was to introduce the parameter plane and 
establish its position among tools available to the engineer. The division 
of servo-system theory into various approaches was arbitrary although it 
seemed reasonable. Similarly critique of present approaches and methods was 
pursued as a means of placing the parameter plane and no claim for complete- 
ness or comprehensive analysis is made. 
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CHAPTER II 



Concepts of the Parameter Plane 

The characteristic equation which is obtained by equating the denominator 
of the closed loop transfer function to zero, is probably the most important 
single factor available for determining system performance. In deference to 
convention the characteristic equation will be designated by 

f(s)= ^ = o (1) 

K=o 

where s is the complex variable of the LaPlace transform 

n is the order of the differential equation describing the system 
are real constants 

However such a designation is an over simplification. Consider the follow- 
ing points. 

1. f(s) which is an equation in the complex domain, can properly 
be regarded as two simultaneous real equations. Namely: 

Real Part ^ 

Imaginary Part (S) = O 

2. The complex variable s is in reality two real variables 
since both its real and imaginary parts must be specified 
before s is fixed. Variables commonly used to specify s 
are 

8 =r 6~ -h ju/ = C + ]//— 

For parameter plane formulations the latter specification 
is the more convenient. 
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3. Invariably a servo system has a number of parameters which 
can readily be varied to suit the designer (i*e. system 
gain) and any such change will affect the coefficients 
For parameter plane purposes it is convenient to assume that 
the a^ are functions of two such variables which are independ- 
ently at the designer’s discretion. That is: ^ ® ^ » /3 ) 

(The existance of more than two such variables does not in- 
validate parameter plane concepts any more than the existance 
of more than one variable invalidated root locus concepts. 
Rather, all but two of the variable parameters should be fixed 
prior to any single application of parameter plane technique). 
With these points in mind a more cogent statement of the 
characteristic equation would be 

k 



f f? ( ^ ) 



Real 



k=o 



> /3 > ^ “ Imaginary 



A/ 






2 ) 



Equations (2) are two equations in four unknowns. Classically the coefficients, 
aj^, were assumed known (the designer assigned values to o(,y5 ) and the pro- 
blem was reduced to one of finding roots for the characteristic equation. 

If these roots were unsatisfactory the designer would assign new values to 
o( and ^ and repeat the procedure. This procedure is somewhat illogical in 
that unknowns (variable system parameters) are initially fixed in expecta- 
tion of obtaining knowns (desired root locations) . A better procedure would 
be to solve for unknowns in terms of the knowns. That is, solve the implicit 
equations (2) for explicit expression of o( and jS • 
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( 3 ) 



Direct determination of system parameters in terms of desired s-plane loca- 
tions 2 and UJ^ is now feasible. However, s-plane roots are not generally 
restricted to a single point; rather requirements are such as to place 
them within a certain region. Suppose that a number of s-plane points on 
the boundary of this region are chosen, that the corresponding values of 
system parameters o< and yS are computed using equations (3) and that these 
values are plotted as points in a cartesian coordinate system. The 

curve connecting these points could reasonably be expected to form a boundary 
for the o<-y0 plane region containing values of system parameters which place 
roots of the characteristic equation within the desired s-plane region. 

(A more detailed interpretation of plane curves is reserved for the 

next chapter). Early investigators were interested in this concept pri- 
marily as an aid to stability analysis and therefore restricted their atten- 
tions to an s-plane contour along the jw axis [4,5]. The subsequent work of 
Mitrovic and others generalized procedures to include constant ^ contours 
but were still limited by restricting and ^ to be the last two co- 
efficients of the characteristic equation [3]. More recently Siljak general- 
ized the theory to allow the coefficients of the characteristic equation to 



rapidly become prohibitive. 

Parameter plane (or o<-^ ) curves can be considered from an alternative 



contain any linear combination of 




[7]. This paper extends the 



procedures further to include the appearance of products in the co- 

p P 

efficients. Still further extension to appearance of and terms and 
beyond is feasible but resulting explicit expressions for and 
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viewpoint which is more suitable to mathematical rigor. Here ( c< , ) and 

( 5 » sets of variables which are related by either implicit 

equations (2) or explicit equations (3). The parameter plane (or c<-yS ) curves 
are then the images (or maps) of the chosen s -plane contours. (The reader is 
cautioned that the mapping is not conformal and basic relations of conformal 
mappings are not preserved) . 



Parameter Plane Equations (Linear Coefficient Case ) 

This formulation is an abbreviated version of that advanced by Siljak 
and as such presents only the salient points [7]. Consider the character- 
istic equation (1) where the coefficients a^ are linear functions of the 
parameters o( and ^ . 

, /8 , s ) = ^ c,( + e^) s =0 

K=0 



( 4 ) 



In this equation Cj^ and are constants and s is the complex variable 
of the LaPlace transform. It is convenient to express s as 



s = LJ. 






(5) 



and 

S* = w; (t, 5) -I- 3 (- KSO W 

where w is the undamped natural frequency and ^ is the damping constant 
n 

which assumes positive values in the second quadrant from zero on the 
imaginary axis to plus one on the negative real axis. Observe that equations 
(5) and (6) are not inconsistent; rather equation (5) serves to define 
and Uj^(- ^ ) as explicit functions of ^ . 

By manipulation of equations (5) and (6) some subsequently useful re- 
lations between the and Tj^'s can be obtained. 
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Argument relations 



\(- 5) = (-1)" \(5 ) (7) 

\(- ^ ) = (-1)’^^' \(^) 

Recurrence relations 

) -2^ V S ) + \.1< ° (8) 

Vi< ? > ^ + \-i< S > ■ ° 



Interrelation 



V 5) - 5 V 



(9) 



Also the Tj^( 2 found to be Chebyshev functions of the 

first and second kind respectively. 

Armed with these relations we proceed to find explicit solutions for 
oC and p . Substitution of definition (6) into the characteristic equa- 
tion as expressed in (4) yields: 

r/ 



The real and imaginary parts must independently equal zero to satisfy this 
equation. 




2 



K*0 



(ocb^ +■ 6,^) UJ„\ (S) 

fic^+ejub'' /r?" uj-s) 



=0 

=0 



Removing common factors and employing relations (7) to change the arguments 
of the Uj^'s and to positive zeta results in the expressions: 

^("0 =0 



^(-0 



/3c,+e^) uC U,(^Uo 
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Equations (9) which express the interrelation between the and can 



be used to eliminate T, . 

k 



2 ^ e.) <J.‘ u , (■ s) - 2 0 * («< /« C K e. ) ^ S) 

JCaO J LK*® 

+ eO wj' U, (5)1 = 0 

_>c»o J 

But the quantity in the left hand bracket is identical to the quantity below 

it which is equal to zero. Incorporating this observation allows us to write 

n 



T (-/)*' (.(b^+zsc.-t-eOH* Uk.,( 5 ) = o 
f’C-')' (c(^+/3c, + eOwJ‘UK(£) = o 



( 10 ) 



Equations (10) are linear in o( and ^ and are easily solved for their explicit 
expression. This point is readily apparent when equations (10) are rewritten 
in more concise notation. 



o( B, ( ^ ) ^v) ^ C| ( ^ ) ^i(^, — O 

where by definition ^ 

B (5 , w„) = y Un.|(5) 

K»0 

B2( ^ ) = J (-/) Uf^ (^) 3 

KsO 

Ci< 5 . U|J., (5)9 ■ H 

CjC ^ ) - ^ (~'fc^ U„ (5) » c 

E 



(11) 



— uJ 



( 12 ) 



K»0 



1<5 .o„) - 

KsO 
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In the interest of brevity the functional dependence of B, C and E will be 
omitted when such omission can cause no confusion. The explicit solutions 
for o( and p obtained from equations (11) are: 

oU^ ,(J^) ^ (13) 

B.C£-6jC, 

/3 ( ^ , U4, ) = Bs E, - 6( E‘g 

B| Cg - QgC, 

These are the explicit expressions for ©(and |S which allow mappings of 

points or contours (except for points on the real axis) from the s-plane 

into the parameter (or o(,-y(? ) plane. Applicability does not extend to 

points on the real axis because the derivation included division by a factor 

k 

of (1 - ^ )^ which is identically zero thereon. 

Relevant parameter plane equations for points on the real axis are 
easily obtained. Define values of the complex variable s on the real axis 
by 

s 4 - (S~ (14) 



For values of s on the real axis the characteristic equation (4) becomes 

K 



f( c< . ^ . <T) = ^ 

K=0 



(15) 



Of particular importance is the degeneration of the characteristic equation 
on the real axis into a single real equation whose imaginary part is identi- 
cally zero. For easier discussion define: 

N 



Br((T) 






K=0 

V 






C|^ (T 



(15a) 



Ej^( (T ) “ 



k=o 
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and write equation (15) as 



c<B„(ir) + ^CR(r) -t-r^Cir) = 0 (i6) 

Equation (16) is a single equation in three unknowns o( i jS and 5" which on 
selection of a real point ) to be mapped defines a curve in the 

plane. For the case under discussion o( and appear linearly and the 
curve will be a straight line. 

Equations (13) and (16) are the keys which allow interpretation of the 
characteristic equation in the parameter plane. Equation (13) is applicable 
for all s-plane points except those on the real axis and it maps points into 
points and s-plane contours into o (-|8 plane contours. Equation (16) represents 
the degenerate case when s assumes a value on the real axis and it maps a 
single s-plane point into a straight line in the plane. The significance 
and interpretation of these parameter curves is explained in detail in the 
next chapter. 

Parameter Plane Equations ( Product Case ) 

All too often system parameters appear in the characteristic equation 
non-linearly; an occurence which precludes parameter plane analysis by the 
techniques thus far developed. Theoretically parameter plane techniques can 
be extended to cover any eventuality but the impracticability if not the im- 
possibility of obtaining explicit solutions for system parameters in terms 
of s-plane root locations prevents theory from becoming practice. Character- 
istic equations in which products appear are a frequently occuring and 

important exception to this statement. 

f(o< , s ) + = o ( 17 ) 

K=0 



13 




■ 4 *^’ 








Parameter plane notation and definitions for the product will be consist- 

ent with those employed for the linear case. Noteworthy is the continuing 
relevance of relations (5) through (9). Separation of equation (17) into 



its real and imaginary parts yields: 

V (-/) -I- /S C|( + Tk C 5) = o 

^ (— /) ^ -I- L(k ~ ® 

IteO 

The same manipulations employed in the linear case allow these equations to 
be expressed as: 

X (-/) ' (o< 4- /s C, + 4- ej U ( 2) = o 



K-O 

N 



(-/)" (x e J a. (5) = o 

K-O ^ 

In the Interests of brevity define: 

B ( 5 . ^1^-' ^ 

K*0 
V 



B ( 5 , w, ) - 2] Ir'f bn U|^(5) 

K-O 

Ci( S . <^1.) = 2 ^K-t 

KsC 

c2(5.w„)-7(-0"c,w;U,(S) 

Di( $ B H J (-O'^ d.^ U,t_, (^) 

D ( ^ , ) = y u{j^ (^) 

^ feo 

VS.H)- $;(-0'eKH'^U,.,(S) 

E2<$ » H ^ Uk (^) 



K-O 



(18) 



(19) 



and rewrite equations (18) as: 

o< 6, ( J , H) C , ( 5 , ) ^- ®<)2 D, ( s , H) f , W,) = o 

0( 6e(5, '*'»)+ ^Cg('5,w„) + + fgCS,H)=0 



(20) 
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(Notice that definitions (19) are identical to definitions (12) with Dj^( 2 > 
) and D^( 2 » ) added) . Equations (20) can be solved for o( and jS 

as functions of ^ and LtJ^ by three different methods. 

1. The parameter is initially eliminated from equations (20) 
and o( is obtained from the single resulting equation using 
the quadratic formula. Values for ^ are subsequently found 
from a linear formula utilizing the previously found values 
of 0< as knowns. Starting with equations (20) 

o<6,-e, 

o(^(B^D2 - B2Dj) +o<(B^C2 + D2E^ - B2C^ - D^E2) + (C2E^ - C^E2) “ 0 



= ( %Q L.t,Q| ^5 ~ S| Cg -Ofi r, ) 
^ e(B,Dz-BzD,), 



( 21 ) 



- teeC. + D,Ep-B.C^-DsF,f- H (8,05-83D>)(C^E, -C,E^) 



e(B,Ds~BsD,) 

' ^ c,^D, + C, 

The last two formulas yield identical results for ^ . 

2. This method duplicates method 1, except that 0( is initially 



/ o( "h C 2 



( 22 ) 



eliminated and jQ is solved via the quadratic formula. 



&(C,LuJ) = (S| Ce -k D| Eg - Sa Ci -Ds£j) 1 

e ( C . D 2 - CgO .) ^ 



(23) 



-UBaC, + DsE:,-BX^-D,Esf - H (C^D^-C^D,)(B^E,-B,E^) 

B(C,Ds-CsD,) 



OJ) = ^ 



(24) 



o<:(2)K,)= - 

. <SD£-i-6a 

3# A third procedure involves solution for both o<( ^ ) and 



^ ^ ) by the quadratic formula after ^ and 0(, 
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respectively are eliminated among equations (20). The re- 
sults are solutions (21) and (23) which will be designated 



«C( $ . UJ^) = XifT 
( ^ ) = a±{T 



While equations (21) through (24) display explicit solutions for ©C 
and ^ a relevant question still remains - are solutions (21) and (22) iden- 
tical to those obtained in (23) and (24)? This question paraphrased to 
apply to the third procedure would be - are each of the possible point sets 
(Oi, t jS ) valid solutions for equations (20) or are two of these sets ex- 
traneous? Direct substitution reveals that the solutions obtained using 
equations (21) and (22) are identical to the solutions obtained using equa- 
tions (23) and (24). Moreover, using the third procedure yields valid solu- 



( o< , ) = (x + fy, u - fv) and ( 0< , ^ ) * (x - \fY, u + fv) 



The two remaining possibilities are extraneous and will not satisfy equations 



As for the linear case s-plane points on the real axis are degenerate; 
the Imaginairy restriction on the characteristic equation is identically 
satisfied and only one equation relates the three variables o(, iS and s. 
This equation which is appropriate for points on the real axis is obtained 
by defining s as in (14) and substituting into characteristic equation 



tions: 



( 20 ). 



(17). 




(25) 



As before define: 




(25a) 




I 



I 



« I 



I Mr 



vr> = S 

and write equation (25^ as: 

/S (o<) (T) = _ o<Bff (<r) -t- Er((T) (26) 

Dff (6“) Cr (^) 

The locus of this equation for any given value of (S~ describes a hyperbola 
in the parameter plane with asymptotes of: 

0< = - ) and /3 = - 

O/iOr) 0 ^((t) 

Comparing these explicit expressions with those occuring for the case 
of coefficients linear ino(^and ^ we find three mentionable differences. 

1. In the linear case all s-plane values mapped into real values 

in the o(-|8 plane. For the product case this is no 

longer true; the possibility for complex values exists. 

Physically this implies that it is impossible to place s- 
plane roots at such locations by adjustment of oC and 

2. In the linear case each complex s-plane point mapped into a 
single point in the parameter plane. For o<^ products a 
complex s-plane point usually maps into two different points 
on the plane. 

3. For both linear and product cases real s-plane points have 
curves as plane images; however, where the curve is a 
straight line for the linear case it becomes a hyperbola for 
the product case. 

Parameter Plane Methods Extended Beyond the Product Case 

Extension of parameter plane methods beyond the product case appears 
prohibitive even when explicit solutions for and ^ can be found. 



17 



Suppose we allow characteristic equations of the form 



N 

f(o< , ^ . S ) + =C> 



it-o 



Under these assumptions explicit solutions for o( and ^ exist and involve 
no insurmountable or extraordinary difficulties; in fact simultaneous solu- 



tion of two 


conics for and ^ 


yields the 


desired expressions. 


o^^Aii + 


0<p ^21 ^ ^31 




^ *51 *61 ' “ 


Aj2 + 


A^2 + ^^^32 


C<A^2 + 


j3*52+ Ajj-O 


where 








^i<£ 


’ ^ ^ ^ 


(5) 



*32< 5 > *"« > Uk (5) 



K=0 

The problem of consequence is interpretation of the parameter plane curves 
since each s-plane point generally has four.o4~^ plane images* Additionally, 
the four valid solutions for ( , ^ ) must be distinguished from among a 
minimum of eight generated in the process of solving the simultaneous equa- 
tions* 



Parameter Plane Methods for More Than Two Parameters 

The parameter plane is a competent two variable technique in analogy to 
the root locus being a capable single variable technique. Pursuing this 
analogy we find that parameter plane technique can accommodate three variables 
with a facility about equal to that of the root locus to accommodate two 
variables. Problems having a third parameter plane variable can be approach- 
ed in three ways. 

1. An additional restriction or specification is invoked (such 
as a steady state error restriction) which permits solution 
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for the third variable, ^ , as a function of the remaining 
two, f( c< , ). With elimination of the third parameter the 
problem reverts to the standard two parameter problem and 
analysis proceeds as before. However linearity of oC and ^ 
in the coefficients of the characteristic equation is 
commonly destroyed when eliminating in favor of f( 0< , 

/3 ) . This shortcoming seriously diminishes this otherwise 
attractive approach since loss of linearity usually precludes 
parameter plane analysis. 

2. Addition of a third dimension to the two dimensional para- 

meter plane allows adequate representation of the three 
variable problem. Under this representation an s-plane 
contour maps into a surface in ( oC » ^ ^ ) parameter 

space. Any choice of system parameters which lie on this 
surface will have at least one pair of roots on the chosen 
s-plane contour. Only the problems innate to three dimensional 
presentation limit the attractiveness of this method. 

3. One selected s-plane contour can be mapped into several o<.~^ 

plane curves, each curve representing a different value of )( , 
the third parameter. This procedure effectively extends para- 
meter plane analysis to three variables especially \dien only 
one or two s-plane contours are of interest. Another credit 
is that applicability is not limited by the manner in which 
the third parameter appears in the characteristic equation. 
Moreover, the plane 'image of any specified s-plane point 

(real axis points excluded) plots on a straight line as the 
third parameter is varied whenever all parameters are linear 
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with respect to one another. 



This last observation is quite useful in the following context. Suppose we 
specify any s-plane point as s(x,y) where x is held constant and y is the 
running variable along the s-plane contour of Interest. The map of this 
contour in the c<-^ plane depends on the third parameter; a different mapping 



represent the same value of the running parameter y - the locus of these 
points is a straight line in the parameter plane. Furthermore, the spacing 



The apparent utility of these linearity properties warrants the more 
rigorous foundation supplied by formal theorem and proof. 

Theorem I . Suppose three variable parameters appear linearly in the 
coefficients of the characteristic equation. 



Then the locus of any fixed s-plane point in the parameter ( ai — ^ ) plane 
is a straight line and the image of the point is a linear function of the 
third parameter 

Proof . Resolve equation (27) into its real and imaginary parts, equate 
each part to zero and simplify as was done when originally formulating the 
parameter plane equations. The result is two equations analogous to equations 



for each different value of ^ . Connect the points which, on each curve. 



of the various image curves is proportional to the value of ^ . 




= 0 (27) 



( 11 ). 



oCB^( ^ ) + /QC^( $ ) + )( F^( 5 , U;^) + E^(^ , ) = 0 



(28) 



0CB2( $ , tu'w ) + /3B2( 5 , 07^ ) + YF2( $ , 07^ ) + E2( ^ = 0 



where 




(28a) 



K = 0 
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K=rO 

and all other quantities are as defined in (12) 

Explicit solution of equations (28) yields the customary parameter plane equa- 
tions. 



rs/ t Y ^ tjJ \ m, y (Ci i-a — CeFi)+ (CiEa— CaEt) 
<x.( y , e , u/„) (b.c^ -BjC,) 

/8or,5,ay«)= y(B^F,-B,F^)■^■(BgF , -e,F^)„ 

— Dg (-i; 



(29) 



For a fixed point in the s-plane ^ and become knowns and oC and ^be- 
come functions of the single variable Y . Under these conditions equations 
(29) are recognizable as the parametric representation of a straight line in 
the oC—jS plane. Furthermore the linear appearance of Y ii' both of these 
explicit expressions implies that the position of the image point on the 
straight line locus is linearly proportional to the value of Y . Q.E.D. 



Parameter Plane Equations in the Z-Plane 

The aforementioned analogy between parameter plane methods in the s and 
z-planes can be strengthened to practical identity of procedures and equa- 
tions. (The only noteworthy difference between the two planes is choice of 
contours to be mapped onto the o<-^ plane; constant $ contours are popular 
in the s-plane whereas spirals of constant damping or circles of constant 
settling time are popular in the z-plane). This close correspondence made 
formulation of the parameter plane equations in the more familiar s-plane 
with subsequent extension to the z-plane the most attractive alternative. 

For z-plane formulation we assume a characteristic equation of the form 

hi 

f(oc, z) = ^ =0 

K=o 
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with provision for inclusion of terms for the product case, three 
parameter case, etc. Results derived in the s-plane are valid in the z-plane 
subject to the following modifications. 

a) If the complex variable z is defined as 

With z so defined no modification of s-plane formulas is 
necessary beyond replacing ^ by and by . 

b) If the complex variable z is defined as 

(?. + W'-Sf ) 

The total effect of so defining z is to obviate the need for 
1/ 

the (-1) factor in the definitions of B,C,D,E and F. 
Specifically, all s-plane formulas are valid if the definitions 
in (12), (15a), (19) (25a) and (28) are redefined as 

A/ ^ 

BjCS* -2. 

K=0 

ic=o 

K=o 

.M,,) = J 

®1^ $2 ® j dj; U-K-l (^Z) 

K=0 

h< U,., ($,) (30) 

K=o 

k=o 



22 



! 





# 



I 






K^o 



) - ^ t, «; 



,w..)=J ?.<U,_,(5.) 

k=o 

-2 

>-Sc« 

K=0 

.-2 j, 

R=o 

>=J 

IC=0 

>-2«."-. 

MO 



<5; 






The reader is cautioned that the analogy between ^ and and and 
extends only to the mathematical formulation of the parameter plane equations. 
The value of is in no way indicative of system overshoot nor is 
indicative of the system's undamped natural frequency. 



Mitrovic's Method and Extensions of Mitrovic's Method 

Although Mitrovic methods are special cases of the parameter plane they 

t f 

warrant mention because of a relation between the ^ (M-plane) curves and 
the s-plane roots not present in the more general case. This relation and 
its application to the z»plane will be discussed in Chapter 3; the present 
objective is simply to define Mitrovic methods. Parameter plane theory assumes 
a characteristic equation of the form; 



f( oC, S ) = O 

| C =:0 

Mitrovic *s method and its extensions assume that only one and one c 
(i ^ j) have non-zero value; oC and ^ are subsequently redefined as: 



j 
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f3 = /3 Cj 

The characteristic equation in a format apropos to Mitrovic formulation then 



appears as: 



If i - 0, j 



f( OC 






s) 



N 



G IX S + 






4 - /3 =0 



K=0 

j 

1 the fomulation is the one originally advocated by Mitrovic; 



otherwise 



the formulation is referred to as an extension of Mitrovic 's method. 
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CHAPTER III 



The Interpretation Problem 

To regard a functional relation between two complex variables, w = g(z), 
as limited to defining w in terms of z is indeed tunnel vision. Lost would 
be such concepts as images of z-plane contours and mappings of z-plane areas 
into selected areas of the w-plane. An equally grievous error would be to 
regard the explicit expressions for o^( ^ , UJ^ ) and ) QS formulas 

to place roots at specified s-plane (z-plane) locations. Such regard would 
obscure two most important interpretations, that of mapping and that of parti- 
tion of the oC~ ^ plane into regions having a similar root property. The 

I 

problem, therefore, is not for a single interpretation of o6-^ plane curves 
but rather it is for an answer to the question, "How can the significant in- 
formation buried in s-plane (z-plane) pole zero configurations be extracted 
from parameter plane curves?" 

Interpretation in the Plane 

A simple yet comprehensive answer to the interpretation problem is pro- 
vided by examining various point sets in the plane for common properties. 
In figure (la) a z-plane contour terminating at two points on the real z- 
axis is designated; figure (lb) shows the oC— ^ plane image of this contour 
and the images of the two real z-plane points terminating this contour. With 
reference to figure (lb) define the following point sets. 

Point Set A - Point set A contains the images of all points on the 
designated z-plane contour. 

Note: This set contains all points on the image curve 

drawn in figure (lb) and no points not on this curve. 
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Point Set B - Point set B contains all points which are the image of 



the real z-plane points 6^^ and 6^^ . 

Point Set C - Point set C contains all points in the cross-hatched area 
and it is bounded by but disjoint from set A. 

Point Set D - Point set D contains all points in the dashed shaded area 
and it is disjoint from set A and set B. 



z-plane 





Let M( oC j ^ ) be any point in the parameter plane and consider the signi- 
ficance of M belonging to certain point sets. 

Observation: If M belongs to set A the characteristic equation must 

have one complex root lying on the chosen z-plane contour and another root 
lying on its reflection in the real axis. Recall that explicit equations 

■/« > 
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give the parameter values needed to place a root at ( ) . 

However, these explicit equations can also be considered as mapping functions 
which establish a relation between z-plane points and points. Regarded 

thus, every point in set A has its domain on the designated z-plane contour 
and satisfies these explicit expressions; hence the conclusion observed is 
inevitable. 

Observation: If M is contained in set B the characteristic equation must 

have at least one real root at z = 0^ (or ) . 

Observation: All M-points located in set C have an equal number of 

characteristic equation roots enclosed by the z-plane contour and its reflec- 
tion in the real axis. Furthermore the number of characteristic equation 
roots enclosed by an M-point in set C differs by two from the number of roots 
enclosed by an M-point in set D. The validity of this statement can be 
established by contradiction. Assume two points exist in set C which have a 
different number of characteristic equation roots enclosed by the z-plane 
contour. Then, by continuous motion of the M-point from one of these points 
to the other, characteristic equation roots could move across the closed z- 
plane contour without ever once having the roots lie on the contour itself. 

(All M-points having roots on the contour belong to set A which is disjoint 
from set C.) This is impossible; therefore by contradiction all points in C 
have the same number of characteristic equation roots enclosed by the z-plane 
contour. Similar continuity arguments can be applied to show that a differ- 
ent number of roots are enclosed for M-points in set C and set D. (To deny 
this proposition would imply that a root lying on the specified z-contour could 
be moved in any direction in the z-plane without affecting the number of roots 
enclosed by the contour). That the enclosed roots in set C and D differ by 
two stems from the restriction that coefficients of the characteristic equa- 
tion be real; root motion across the contour’s upper limb requires motion of 
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the conjugate root across the contour's lower limb. 

Observation: The plane images of points at which the z-plane 

contour intersects the real z-axis also serve to divide the parameter plane 
into point sets. M-points belonging to sets on opposite sides of such con- 
tours will differ by one in the number of characteristic equation roots en- 
closed by the z-plane contour. Justification of this statement is analogous 
to that of the preceding observation. 

Preceding observations although they cannot claim mathematical rigor 
provide a framework for parameter plane interpretation one step above the 
intuitive approach usually advanced. For application these observations 
can be profitably restated. 

1. M-points on an image curve imply a characteristic equation 
root on the z-plane (s -plane) contour upon whose image the 
M-point lies 

2. Parameter plane images of z-plane contours and images of 
the intersection points of these contours with the real 
axis serve to divide the ^ plane into regions or point 
sets. The distinguishing characteristic of each region is 
that all M-points within a region have the same mamber of 
characteristic equation roots enclosed by the z-plane (s- 
plane) contour. 

3. If the boundary between regions is the image of a complex 
contour the number of roots enclosed for M-points located 

in adjacent regions differs by two. If the boundary between 
regions is the image of a real axis point the number of roots 
enclosed for M-points in adjacent regions differs by one. 
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Shading Convention and Rules 

The interpretation of curves as images of selected s-plane (z- 

plane) countours was previously discussed. Moreover, the preceding para- 
graph showed that movement of the M-point across an curve corres- 

ponds to two complex roots moving across the selected s-plane (z-plane) 
contour. The direction in which the roots move across the s plane contour 
is usually available by inspection of the direction in which the M-point 
moves across the curve. (i.e., Motion of the M-point across the o<-/5 

plane image of the contour ^ = 5^ in the direction of oi—jQ plane images of 
contours $ implies that the s-plane roots move across the contour in 

the direction of increasing ^). However, parameter plane regions in which 
images of two different ^ contours intersect or in which the image of a 
given s-plane contour loops on itself do not readily allow the direction 
of root motion across the s-plane contour to be inferred from the motion of 
the M-point across its o4-y8 plane image. For such cases, and in fact when- 
ever required, the direction of root motion across an s-plane contour as 
the M-point moves across its image in the parameter plane can be deter- 
mined mathematically. 

Adopt the convention that the parameter plane images of constant ^ 

(or constant ) s-plane contours are shaded so that an M-point cross- 
ing the o6-^ curve from shaded to unshaded side causes two complex s-plane 

roots to leave that portion of the s-plane enclosed by the contour. 

Theorem I The sign of the Jacobian J* f ) determines the side 

of the image curve to be shaded. If the Jacobian is positive (negative) 
the left (right) hand side of the image curve should be shaded facing in 
the direction of increasing i ^ )• 

Proof Construct a cartesian three dimensional vector space on the 
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parameter plane whose position vector is: 



r(oC,/<3.Z) = oCx-h/3j+2^ 

In this space all curves lying in the parameter plane, as all images of s- 
plane contours must, will have position vectors of the form 

r(oC,^)= 

Now recall that a requirement for the construction of parameter plane curves 
was the existence of explicit solutions for oC and |3 in terms of ^ and 
UJ^ . Namely: 

C)C - ©<.(?* ) 

/S. 

This requirement allows the position vector oC, , ^ ) to be written as 

r(f,uj^)= oC( $ , ^ 

The oC-yS plane image of a constant ^ contour of value is traced by 
the vector 



- oC(p^, /3( 



Fix attention at one point on this curve and determine the motion of that 
point resulting from an incremental increase in ^ .It is 

o)^ 

Obviously the direction of M-point motion in the direction of increasing ^ 

at this point is determined by the vector Since an M-point 

/• 

moving in a direction to increase > causes roots to enter the area en- 
closed by a constant $ contour that side of the o^— plane image toward 
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which the vector points is shaded. Graphically this scheme 

suffices but it can also be stated in a mathematical formulation. By 

definition the tangent to the plane image of the constant s- 

plane contour at the point ( , IaJq ) in the direction of increasing 

UJ^ is given by the vector . If the counter-clockwise rotation 

of the vector into the vector is more than 180® the 

c)uJ^ 

latter vector points to the right of the constant curve image and that 
side of the curve should be shaded. Contrariwise, if the rotation is 

less than 180® the vector ^ points to the left of the constant 

mapping in the o4— ^ plane and that side of the curve should be shaded. 
Fortunately this criterion is the exact definition of the cross-product of 
two vectors lying in a cartesian plane. 






If A is negative shade the right hand side of the plane curve facing 

in the direction of increasing \jJ^ ; if it is positive shade the left hand 
side. At this point reconcile the sign of A with a more familiar expression. 
To do this express r( ^ , UJ^ ) in terms of its oC and ^ components and 

take the required derivatives* 



c^r ( Co, UJo) _ c^oC ^ j ^ 

<^UJ^ 

c9r { ^o, ^o) _ 1 ^ I . UJe) A 



and ^ 

\/ c)rfa,U;o) 






cPoc ( Co , (*^o) 



o)u^ 



k 
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Hence the sign of A depends on the sign of the determinant which is by 
definition the well known Jacobian 

Extension of these results to other s-plane contours and also to z- 
plane contours is possible. In fact some general theorems must exist 
which encompasses all alternatives, but the statement of such a theorem and 
the determination of necessary and sufficient conditions are difficult. 

Since there are but a limited number of commonly used s-plane (z-plane) 
contours it suffices to list the governing Jacobian for each. The proof in 
each case proceeds exactly as that presented above. 

Extension of Theorem I The sign of the Jacobian T (tf ) , where 
the independent variables x,y are dependent on the contour being mapped, 
determines the side of the image curve in the plane to be shaded. 

Shading rules are as listed in theorem I. 

a) For constant settling time or frequency contours in the s- 



plane the governing Jacobian is 









c)/^ ^ c)/S c)^ 



b) 






where s = (- $ -f- -j ^ | ~ ^ ^UJ 

For constant C or (jJ, (settling time) contours in the z- 

plane the governing Jacobian is 

c9oc ^/s 

where 2 = Ct/ ( - ^ -i- j |// ~ ) 
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Note: Frequently z is expressed as z = {/J^ \^4 

in which case the sign of the Jacobisn is reversed, 
c) For constant damping contours in the z-plane the governing 

<)K^ cPuJ^ 

S C^e ^ ^ 

o^C <^Ci ' 

Mathematical formulation for shading is virtually complete except for 
the degenerate case of points on the real axis. Each of these points have 
an entire curve as their image in the parameter plane since only one equa- 
tion exists (16) which relates the three variables oC, ^ , 5" . (The restric- 
tion that the imaginary part of the characteristic equation equal zero is 
identically satisfied by points on the real axis). If parameter plane curves 
can be constructed explicit solution of equation (16) for oC in terms of ^ 
and (T or for y8 in terms of oC and exists. Designate these explicit 
expressions as 



Xtl)- 



Jacobian is 






















c)oC 

















(oC,(T) 
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Theorem II The shading of the parameter plane curve \rfiich is the 

image of a real point (J'( ) in the s-plane (z-plane) is determined by 

the sign of the partial derivative ^ (or i ) where 

<^T <9 <T 

( ) is a point on the parameter plane image of . The shading 

rule is as follows. 

(a) If Increasing (T ( or 0^ ) from 0^ causes a real root 

to enter the area enclosed by the contour shade the left 
(right) side of the curve facing in the direction of increas- 
ing OC when 

(9(t 
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is positive (negative) . 



Equivalent to this is the statement that if increasing 5“ 
causes a real root to enter the area enclosed by the contour. 



shade the left (right) hand side of the curve facing in the 

direction of increasing fi when is negative 

' o?«r 

(positive) . 

(b) If increasing (S' (or (5^ ) from (S^ causes a real root to 

leave the area enclosed by the s-plane (or z-plane) contour 

shade the right (left) side of the oC-jS curve facing in the 

direction of increasing oC when is positive 

cP(T 

(negative) . 

For comparison a table presents this theorem in a most effective manner. 



Case 


face 

Direcnon 


S/qn of 
qlSK, (S;) 


S^ade 


Face 

Direcr/on 


S/g/i of 

^oC (Jo) 


Shade 


c)o- 


< 9 (T 


Increased (T 
(or 05) implies 
a. rooT enrers 
f/»c. enclosed 
area 


IncreaS- 

mg 

oC 


Pos / five 


Left 


Increas- 

1A9 

/? 


A/egative 


Left 


A/ega^ive 


Ri^hf 


Positive 


Riqff 


Increased (T 
(or 0^) implies 
a, roof leaves 
I'he enciosed 
area 


Increas- 

109 

od 


Positive 


Right 


IncreciS- 

jng 

/s 


Negative 


Plight 


Ne^ahve 


Le^f 


Positive 


Leff 



Table I 
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Values on the real s-axis and the real z-axis are not defined identically; 
the usual case being 

Real s = — (T" 

Real z = 4- 

The manner of definition does not affect the validity of the theorem if the 

user is consistent throughout. For instance, under these definitions a 

change of s from -1 to -2 would indicate an increase in 6^ whereas a 

similar change in z would indicate as a decrease in 0“^ . Notice also 

that either o/Oo) sufficient information to 

cP(T ^<r ^ 

shade and their dictates must always agree. An exception occurs when the 
tangent to the curve is vertical or horizontal in which case one of 

the partials will degenerate to a value of zero or infinity while the other 
will still impart useful information. 

The proof to this theorem is analogous to that presented for theorem I. 

In this case either oC and (T' or j3 and can be considered 

independent variables depending on which variable ( OL or |S ) is 
explicitly expressed. The position vectors are 

P(oC,C^) = odX + explicit expression of 

p ( |3 , ) = oC. + explicit expression of oL 

From this point forward the proof is identical with that of theorem I. 

A few concluding comments on shading and shading formulas should be 
reiterated for the reader. First, the shading of parameter plane curves 
provides a comprehensive picture of root motion across selected s-plane 
contours as a function of the simultaneous variation of two system para- 
meters. When several s-plane contours are mapped onto the plane this 

root motion is usually evident without the tedium of shading and if shading 
is desired it can be accomplished by inspection. For any case where inspec- 
tion shading is impossible the formulas presented here resolve the difficulty. 
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In this they extend beyond existing criteria which can be applied only for 
coefficients a^ linear in oC and ^ [4]. However, the labor involved 

in evaluating the sign of a Jacobian is somewhat greater than that required 
for existing shading procedures and the former undoubtedly will find employ- 
ment only as a last resort. Probably the most significant contribution of 
these rules is that they put shading on a more rigorous framework. 

Example of Parameter Plane Curve Interpretation . 

An example will serve to illustrate the salient features of parameter 
plane interpretation and curve shading discussed in the preceding sections. 
Consider the third-order error-sampled system with zero-order-hold and tacho- 
meter feedback shown in figure 2. Suppose that there is a requirement that 
roots for the sampled transfer function have zeta values greater than 0.5 or 
lie within the corresponding spiral contour in the z-plane. Allowable s- 
plane and z-plane root areas are shaded in figures 3a and 3b respectively. 
Proper adjustments of k and k^ to meet this requirement are available from 
interpretation of the parameter plane image of the z-plane contour which 
bounds the allowable roots. Choose oL and ^ as k and kk^ respective- 
ly; for this choice the characteristic equation of the system shown in figure 
2 is: 

f( ^ , z) = z^ + (.058 oC + ,310/G - 1.974)z^ + (32) 

(.16304 - .122 /S + 1.198)z + (.028o4 - .188 /S -.223) = O 

The coefficients of this characteristic equation are linear in OC and ^ ; 

plane curves (mappings) are obtained by utilizing the parameter plane 

formulas for the linear coefficient case developed in chapter 2. These 

image curves and the images of the two points at which the designated z- 

plane contour intersects the real axis are shown in figure 4. Shading in 

this figure is by inspection. (In this case inspection shading can be 
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R(s) 




Figure 2 




M- Point 


Values 

i 

Alpha Beta 


Characteristic Equation 
Roots 


Number of 
Roots Within 
Contour 


Ml 






.5I1-2 


± j369 


.712 


3 


M2 




.1 


.808 


t j326 


.30i+ 


1 


M3 


-.1 




.389 


± j392 


1.0^7 


2 


Mlt- 


-.1 


.75 


.355 


± jll-77 


1-039 


0 




.3 


.5 


.1+97 


± j368 


.808 


3 


m6 


.5 


.7 


.I1.96 


± ji}-66 


.736 


1 
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achieved as follows. Recognize that the characteristic equation is of order 
three in z; further recognize that for ^ ^ 0 the roots to the character- 
istic equation are identical to the open loop poles which are known to be 
within the contour. Hence all three roots to the characteristic equation lie 
within the designated contour for M-points located in the region which in- 
cludes the parameter plane origin. Since this is the maximum number of roots 
which can be enclosed an M-point traveling across any boundary to this region 
implies roots leaving the chosen z-plane contour. With these observations 
shading can be effected without recourse to formulas). Once the para- 
meter plane curves are drawn and shaded (as in figure 4) the compensation 
problem reduces to locating the M-point in the region which results in all 
roots lying within the specified z-plane contour. The preceding discussion 
indicates that this region is the one which includes the parameter plane 
origin. To verify this analysis six M-points in four different enclosed 
root regions have been indicated in figure 4^ and the roots to the character 
istic equation for each choice have been computed. Table 2 lists these roots 
and indicates how many fall within the specified z-plane contour. For each 
M-point theory and practice agree. 

Let it be emphasized that this example illustrates parameter plane inter 
pretation and shading procedures rather than servo-system design in the para 
meter plane. For the latter more contours (images) are drawn; shading may 
or may not be required and additional information to guide selection of an 
M-point may be displayed. 

Mitrovic Curve Interpretation for Common Z-Flane Contours 

Mitrovic methods (defined in chapter 2) are special cases of the para- 
meter plane which warrant attention in linear (s-plane) theory because the 
Cauchy criterion regarding argument change for roots enclosed by s-plane 
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contours is directly interpretable as encirclements of the M-point. For 
the Mitrovic plane the tedious sub-division of the parameter plane into accept- 
able and inacceptable regions can be bypassed in favor of a simple inspection 
which determines whether the Cauchy criterion for enclosed roots is or is not 
complied with. The relation between constant zeta curves on the Mitrovic 
plane and the Cauchy criterion has been thoroughly explored in the literature. 
(6, 8). A compendium of these relations is simply this. When all roots of 
the characteristic equation are enclosed by the s-plane contour being mapped, 
argument f(s) where f(s) becomes the characteristic equation when equated to 
zero, increases monotonically . In turn this requires that the M-point be en- 
circled in a predictable manner by the Mitrovic plane image of the s-plane 
contour. Were constant zeta contours meaningful in the z-plane, s-plane 
interpretation would be applicable to them; however, interesting z-plane 
contours are spirals of constant damping and circles of constant settling 
time. It can be shown that when all roots of the characteristic equation 
are within either of these contours argument f(z) increases monotonically 
as the contour is traversed in a positive direction. (Proof can be obtained 
by resolving f(z) into its factors and showing that each has a monotonically 
increasing argument if it is located within the contour). At this point 
complete analogy between s-plane and z-plane interpretation of Mitrovic 
curves fails. The monotonic increase in argument still sllows satisfaction 
of the Cauchy criterion to be determined from the manner in which image curves 
encircle the M-point in the Mitrovic plane. However, a few simple 

examples suffice to show that this interpretation does not extend to general- 
ized Mitrovic coefficients for z-plane contours of constant damping and con- 
stant settling time. 
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Recall that the characteristic equation and Mitrovic polynomial are 



defined as: 



f(z) = F(z) + aj^z + 



f^(z) = F(z) + B^(z)z + Bq(z) = 0 



Subtraction of zero (the Mitrovic polynomial) from the first equation yields 



f(z) =• (a^ - B^(z))z + (a^ - B^(z)) 

Consider the f(z) locus and its Mitrovic ®q) curve shown in figures 5a 



and 5b respectively. 




BqCz) 



Re[f(z)l 




Bj^(z) 



Figure 5a 

If all zeros of f(z) are enclosed by the z-plane contour the f(z) locus al- 
ways starts on the positive real axis and argument f(z) increases monotoni- 
cally as the locus f(z) encircles the origin. The total change in argument 
f(z) is 7T radians times the number of roots of f(z) enclosed by the z- 
plane contour. Since argument f(z) is a continuous function (assuming no 
root of f(z) lies on the specified z-plane contour) the Cauchy criterion re- 
quires that 



argument f(z^) = kT k = 1, 2, 3, , n-1 



41 






■4 




I 



* • 



(that is f (z^) be real) for n-1 values of z on the specified contour where 
f(z) is a polynomial of degree n. the n-1 values of z do not include the 
two points at which the specified contour intersects the real axis. At 
these two points f(z) must also be real but the interpretation which follows 
will not be valid. Take particular note that under these restrictions the 
imaginary part of z^ is not zero at these n-1 points. Hence 

f(z) = - B^(z)] z^ + [Hq - Bq(z^)] 

being real requires that real. This occurs only if a^^- 

Bj^(z^) = 0; which in the Mitrovic plane requires that the curve lie 

on the vertical line through the M-point. Furthermore since argument f(z) 
is monotonically increasing the half lines aQ">BQ and a^^^B^ must be cut 
alternately with the half line a^ < B^ cut first. That is: 

argximent f(z) = 2kTf implies a^ >Bq 

argument f(z) = (2k + 1)7T implies a^ <Bq 

The half line a^ <Bq must be cut first because argument f(z) is monotonic 
and must equal TT , a point on the half line a^^BQ, before it can equal 
2TT , the point at which it crosses the half line a^^BQ. 

The points at which the specified z-plane contour intersects the real 
axis were excluded because previous arguments did not apply to them. Observe 
that if z^ is real 

f(z.) = [a^ -B^(zp] z. +[ap -Bq(z.)] i* real 

and there is no requirement that ~ general this term is 

not equal to zero for real z and this fact can be used to determine the manner 
in which the B^ -B, curve encircles the M-point. If z is the value at which 
the specified z-plane contour intersects the positive real axis then: 
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«1<V > “I 

“i<v < =1 



implies counter clockwise encirclement 
implies clockwise encirclement 



This follows directly from the requirement that the line must be 

cut first. The value be computed by the relatively simple formula 



B, 



K=H 



Mitrovic curve interpretation for z-plane contours can be summarized. 

1. If all roots of f(z) are located within a specified z-plane 
contour (settling circles or damping spirals) the B^(z) - 
Bj^(z) locus encircles the M-point in the Mitrovic plane. 



The direction of encirclement is determined by: 

implies counter clockwise encirclement 
implies clockwise encirclement 



Bj(=r)>n 



3. The line a^ - Bq = 0 must be intersected (n-1) times plus 
one additional time each if the Bq - B^ curve starts or ends 
on this line. 

In linear systems the next attempt was to extend Mitrovic 's method to 
include any two coefficients of the characteristic equation. Desirable regions 
in the generalized plane were again found by inspection in the same manner 
as for the B^ - B^^ case. No such simple interpretation exists for z-trans- 
forms when generalized Mitrovic coordinates are used. Since zeta varies in 
all z-plane contours of interest the requirement that the B^^-B^ locus en- 
circle the M-point (a^-a^) is no longer applicable, (i.e., the Cauchy criter- 
ion can no longer be interpreted in the B^-B^ plane by inspection). A few 

illustrations will suffice to illustrate the point. 

2 



Illustration I: f(z) = z + j^a^-B^Cz)^ 
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The criterion "f(z) real" is met when z = j (without restrictions on 



the terms in brackets). Thus alternate cutting of the half lines 
a^ and B2 *^^2 M-point nedd not be encircled 

to meet the requirement that argument f(z) be monotonically increasing. 
Illustration II: f(z) = [ z+ ^aj^-Bj^(z)j ] z^ 



The criterion "f(z) real" is satisfied if argument z equals minus 

argument [ i argument 

requirement on f(z) cannot be determined by simple inspection of the 

B, , , -B, curve . 
k +1 k 

Although it is probable that criteria for desired Mitrovic plane regions for 
generalized coefficients can be found, they will be complicated and it is 
recommended that parameter plane interpretation methods be used. And, with- 
out simple interpretation methods the Mitrovic plane loses its special 
attraction over the parameter plane and reverts to simply a special case of 
the latter. 
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CHAPTER IV 



Preceding chapters endeavored to present and Interpret the character- 
istic equation In a meaningful manner In the parameter plane. Parameter 
plane presentations need not be restricted to the characteristic equation; 
In a more general consideration the parameter plane Is a space In which a 
system Is described In terms of Its two variable parameters. In analogy 
to characteristic equation presentation, system description In the para- 
meter plane Is the reverse of descriptions normally seen. Parameter plane 
presentation displays suitable values of oL and to meet a specified 
system performance criterion whereas conventional procedures first fix oc 
and ^ and display system performance resulting from this choice. Signi- 
ficant advantages acrue when system criteria are described In the parameter 
plane. 

1. The usual "trial and error" search of parameter settings to 
meet specifications Is eliminated. This advantage Is espe- 
cially attractive when specifications are of the form "equal 
to or less (greater) than" and as such are not tractable 
mathematically. In these Instances the locus satisfying 
equality usually divides the parameter plane Into accept- 
able and unacceptable regions where parameter settings can 
be chosen by Inspection to satisfy both this and other 
criteria. 

2. Simultaneous consideration of all criteria and character- 
istic equation root Information Is accomplished by Inspec- 
tion. Tacitly assumed Is that each criterion plots as a 
curve which divides the parameter plane Into regions suit- 
able or unsuitable for that criterion. Parameter settings 
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within the intersection of all suitable regions simultaneous- 
ly meet all specifications. 

Liabilities are also present in parameter plane presentations although they 
might more properly be classified as limitations. 

1. System descriptions are limited to description in two para- 
meters (even though additional parameters may be available 
for adjustment) . 

2. Specified criteria may be difficult or impossible to deter- 

mine as functions of adjustable system parameters. A maxi- 
mum overshoot specification^ typifies such criteria. Theoreti- 
cally a locus of plane points which have the specified 

maximum overshoot exists but in practice the locus is 
virtually impossible to determine. Time-integral criteria 
and magnitude of frequency response at resonance also seem 

to fall within this class. 

These limitations not withstanding, two very important performance specifi- 
cations, steady state error and bandwidth, can be displayed on the para- 
meter plane. They complement parameter plane information already avail- 
able from characteristic equation roots by relating system performance to 
the zeros of the closed loop transfer function as well as to the poles. 

(The formal consideration of closed loop zeros is almost universally neg- 
lected by other transform analysis techniques). 

Steady State Error Curves in the Parameter Plane 

The steady state error or inability of a servo-system to follow an 
aperiodic input disturbance is an important criterion and frequent system 
specification. For linear systems or linear sampled data systems steady 
state error is usually obtained by application of the final value theorem 
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to the LaPlace or z- transform of the error function. This theorem also 
allows the requisite steady state error curves to be computed for para- 
meter plane display. Derivation of the equation governing these curves in 
the parameter plane will be for the z-transformation with a subsequent 
paragraph Indicating modifications required to extend applicability to the 
LaPlace transformation and continuous systems. 

Assume that the open loop transfer function of a unity feedback system 
(or a system which has been reduced to unity feedback form) is 

(z-1)-^ D^<s) 

where D (z) contains no factors of (z-1) nor factors in common with N (z) . 
o o 

Further assume that both N (z) and D (z) contain no terms above products 

o o 

in the parameters oC and . These assumptions are illustrated in figure 
6 and equations (33) . 







K=o 



Do<z) - J (c?c 4- /<3 + oC/3d^ Ck) 

K=0 



(33) 



Typical steady state error specifications require aperiodic inputs of steps, 
ramps, accelerations, whose z-transforms have the form 



R(z) - 



(z-ir 



-1 



where R*(z) contains no factor of (z-1) or (z-1) 

m = 1 Implies step inputs, m = 2 implies ramp inputs and 
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m = 3 implies acceleration inputs. 

The 2 -transformed error signal to which the final value theorem will be 
applied is obtained by block diagram manipulation. 



E(z) 



Do(z)R*^(z) 
>Vo(z)+ {z-\y Do(z) 



The final value theorem then gives the steady state error as 



E 

s • s • 



I/M 

z -»i 




D.(z) R*(z) 

(z-l)"Do(2) 



(34) 



Consideration of equation (34) is divided into three cases. 

1. J?+l-m^0 In this event the factor of (z-1) or (z-1) ^ 

present approaches zero or infinity as the limit is passed and 
error specifications have no meaning for the input applied. 

J? + 1 - m = 0 and = 0 This denotes a type zero system in 



2 . 



which the steady state error for a step input is 



s. s. 



DoO) R 0 ) 

KO)i-Do(i) 



(35) 



3. ^+l-m=0 and 4 0 This denotes type I systems and above 

in which a factor approaching zero as a limit multiplies D^(l) 
in the denominator. The steady state error expression appropri- 
ate here is 

. D. (i)R* 0 ) 

NoO) 



E 



s. s. 



(36) 



Although their present statement does not so indicate equations (35) and 
(36) are the relations which permit steady state error specifications to 
be drawn on the parameter plane. Statements adapted to parameter plane 
use can be obtained by evaluating D^(l) and N^(l) in equations (33) and 
defining a more concise notation for the resulting sums. That is: 
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1 

1 . 




Nod) 

Do(l) 



^ (o^ /k -H ^ 

MO 

2 (ocb; + /s< H- «/ej; ^- e;) 

k=o 



Now define: 




With these definitions N (1) 

o 

N^(l) = Oi Ej, + ,3 Ep + 

Do(l) = Eg + ^ + 



'p ■ Pk 

KsO 




KrO 



can be expressed as 
®E 



(37) 



(38) 



where all E% are constants. Useable expressions for steady state error 
curves on the oC- j3 plane are rendered by substitution of expression (38) 
into equations (35) and (36). For substitution into (36) 



E s "f" ^ R (0 

■ (otErp+/3E,*»<^r5+E,) 

In this expression all E's are constants as is R*(l). When E is equated 

s « s • 

to the steady state error specification the equation describes the locus of 
all oC-j8 plane points exactly meeting the requirement. The locus explicit 
ly expressed for /3 as a function of oC is 



^ [Ee - Ep E,.J + fEt R"o)- Eh C..i1 
»c[Eo/?"i:0 - Eo,fw.] + [EcR"(i) - 



(39) 



For type I systems and above equation (39) defines the locus of all points 

(OC»/S ) having a steady state error E to an input whose z-transform 
I s • s • 
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is R(S)_ . 

For type zero systems an analogous formula which describes steady state 
error loci for an input step of magnitude R can be obtained by initially 
substituting expressions (38) into equation (35). The final result is 



13 





[Eg ( R~Es.s.) ~ EpEs.^ 


kl 


[e e (R~Es.s^ ~ Eh 


1 




[^£p(f?-Es s.) - ^ Es.sJ 


kl 


[£cO?-e:s.O-EqEs.3 



(40) 



Observe that the steady state error loci given by equations (39) and (40) 
describe hyperbolas in the parameter plane for the product case. 

For the linear case = 0 and the steady state error loci become 

straight lines in the parameter plane. 



Steady State Error Loci for LaPlace Transformed Variables 

If anything steady state error loci are easier to obtain from error 
signals described by their LaPlace transform as opposed to signals des- 
cribed by their z-transform. As for the sampled data case assume a known 
open loop transfer function with numerator and denominator linear in oL 
and jS but not necessarily linear with respect to each other. That is: 

) s (41) 

DqCs) = ^ (c>C -t- ^ 

K=0 



\^ere does not contain s as a factor and the complete denominator 



is 



D 



(s) 



Again consider only inputs of the step, ramp, acceleration type. 



Input = Rt 

Under these assumptions and definitions th^e LaPlace transform of the error 
signal is: 

E(s) = 



S^Do(S) (m-i)IR 



Do(S) 



For stable systems the final value theorem can be applied to exhibit steady 



state error as: 



s.s. 



//Al 

S — 



[sE(s^ = 



/l A\ 

s 



Dq^'S) 

K (s) -h OoCs) 



(42) 



On passing the limit results can again be divided into three cases. 

1. + 1-m 0 In this case steady state error is either zero or 

infinite and no finite adjustment of the parameters can alter 
this result. 

2. ^ + 1-m = 0, 0 These requirements describe a system of type 

I or higher. The parameter plane locus of points having a steady 

state error of E as obtained from equations (41) and (42) 

s. s. 

is: 



04 


[b; («-/)' 


k! 








M 


[c («-/)//? - 



(43) 

3. 1-m =0, = 0 This requirement is indicative of type zero 

systems which are characterized by having steady state errors to 

step inputs. The parameter plane loci for a steady state error 

E to an input step of magnitude R are governed by the equa- 
s . s . 



tion: 



-1 


[L(r-e,,,)-p:e,2 


1 + 






04 I 




k 







(44) 



Observe that application of the final value theorem makes only those terms 
in the numerator and denominator not multiplied by s (after s factors in 
manerator and denominator are removed) significant in determination of 
steady state error. If parameters oC and /6 do not appear in these terms 
steady state error is not influenced by their adjustment regardless of where 
they appear elsewhere. Furthermore, only the manner in which and |S 
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1 ^ 




I 







4 



f 



appear in these terms affects the shape of the steady state error curves. 

If oijS products appear in these terms the steady state error loci will 
be hyperbolas; if o4 and ^ appear linearly the loci will be straight lines. 

Bandwidth Curves in the Parameter Plane ' 

The bandwidth is significant because it indicates rise time or speed 
of response, it measures in part the ability of the system to reproduce 
the input signal, and it approximately describes the filtering character- 
istics of the system.^ 

Bandwidth specifications can be presented on the parameter plane with no 
restrictions above those inherent to explicit solution of oC and yG 
When so displayed, bandwidth can be considered in conjunction with other 
specifications presentable on the parameter plane and satisfaction of all is 
virtually by Inspection, if such satisfaction is possible. Before embarking 
on the derivation of formulas which display bandwidth criteria on the para- 
meter plane some preliminaries deserve attention. Among these is agreement 
on a definition for bandwidth. 

Excite a servo-system by a sinusoidal test signal and form the ratio of 
peak fundamental component output to peak input. (If the system is 
linear the output consists wholly of a fundamental component). The 
lowest frequency at which this ratio falls below a designated value A 
is defined as the bandwidth of the servo-system. 

The most widely accepted value for A is 0.707 (-3 db) although unity and 
other values have been used. A readily accessible expression from which the 
ratio of peak fundamental output to peak input can be computed is the closed 
loop transfer function. If Fourier transforms are employed the magnitude 
of the closed loop transfer function represents this ratio. The Identicdl 
ratio is obtainable from the magnitude of the closed loop transfer function 

^Truxal, J. G. Control Systems Synthesis. McGraw-Hill, 1955: 77. 



52 



expressed in terms of the complex LaPlace variable s if the substitution 
s = jw is effected. For sampled signals a similar ratio is obtained by 
examining the closed loop z-transfer function for values of z on the unit 
circle. 

Parameter plane methods are not primarily concerned with a determina- 
tion of this ratio; rather the values of system variables (locus of 
plane points) which meet bandwidth specifications (cause the output-input 
ratio to be of proper magnitude at the specified frequency) are sought. 
Functionally the parameter plane bandwidth problem is expressed: 






uJbw) 



(45) 



where by adjustment of system parameters oC and the magnitude of the 
closed loop transfer function is forced to have a magnitude A at the band- 
width frequency w^^ ^ ♦ Equation (45) while sufficient to define bandwidth 
loci on the parameter plane is not conveniently utilized for their construc- 
tion, The remainder of the section is devoted to remedy of this problem. 

In line with parameter plane assumptions the closed loop system trans- 
fer function is required to have the form 



L(2> _ K=o 



(46) 



In the course of the derivation it proves convenient to separate both the 
numerator and the denominator into their real and imaginary components. 
The Chebyshev functions introduced in equations (5) through (9) provide a 
satisfactory vehicle for this purpose. For reference, z defined in terms 
of the Chebyshev functions is 

* ~ \_^Z. ^ ^ // ~ 
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To promote brevity the real and imaginary parts of each sum are defined as: 



Rp( 5. . ) - X Pk Tk ( c.) 

V i/'- if Uk(0 

^ K^O 

A/ 

V > “ X (ij) <^7) 

K =0 

»/ 

K =0 



The remaining twelve quantities in this set are similarly defined. It is 
noteworthy that by this definition all quantities R and 1 are real. On 
replacing the sums in equation (46) by their defined equivalents in equa- 
tion (47) one obtains: 



CCt.) _ ix (^p 

R(2) ~ x(Rg4-]Ie)+/J(Rc + 3lc) + <^('?£.^-3lo>(l?E<-3lF) 



Bandwidth specification that the ratio of sampled output to sampled input 
have some magnitude A at a bandwidth frequency of Wj^ ^ equates to the re- 
quirement that C(z, ) over R(z, ) have magnitude A where z_ is the 

D*w« b*w* b«w* 

value of z on the unit circle which images the bandwidth frequency s « jw, 

o • w 

in the s-plane. When these bandwidth conditions are introduced the quanti- 
ties R and I become constants and equation (48) describes the locus of all 
points (oC , ) which satisfy system bandwidth specifications. That is: 



A = 



A 



2 



C(Zt.u,.) 

R (Zi.w.) g 

(oC Rq->- + 

(•^ f?B+ /QRc 



where A is the magnitude of attenuation specified at the bandwidth frequency. 
All R's and l‘s are constants. 



Performing indicated multiplications allows this equation to be render- 
ed as: 
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-I- 2oC^/S(RsR,^IeI^) + 2<^/e^(R^Ro->-Z^lJ + 

o<.^(Rl^ll)+- (RiRc^Rt^R^^IsIc+IoW-h /32(/^2 4- i |) 4_ 



2o^(RpRh •*• Jp Jh) ~*~ (R^Rh* IqIh) (/?H ~*~ Th ) 
eoc(/?e/?,4-X6l^) + e/«(/?c^£+IcTe)+ (Rl-nf) 



Equation (49) describes the locus of parameter plane points which satisfy 
a given bandwidth specification demanding an attenuation A at bandwidth 
frequency w. . Its appearance is forbidding but after recognition of 

D • W • 

certain features the equation itself is not formidable. Foremost is that 
o4 and ^ individually appear to no power higher than two. Thus to find 
a point on this bandwidth locus for a given alpha value Involves no more 
than solution for y6 by the quadratic formula. (Computer plotting of such 
a curve is trivial). It is also significant that equation (49) describes 
a conic section whenever the numerator and denominator coefficients are 
linear in and /S in which case = 1^ = = 0. 

To demonstrate the first of these observations explicit solution for ^ 
will be effected under the assumption that oC is known (assigned). For 
convenience first define: 



Xj( oO » (R^ + I^) + 2 oi (R^Rg + I^Ig) + (r2 + 1^) 

X2(oi ) = (R^ + Io> + 2 OC(R^R^ + Ij,y + (r 2 + l2) 

Yj(OC) = + 2^(RjR^ + R^R^ + + 



( 50 ) 
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) = 2 + I^I^) + 2 C< (R^R^ + Rj^R^ + 1^1^ + l^l^y + 

2<‘>c'^ ^ W 

W^(0<. ) = Od (R^ + I^) + 2 OC(RpRjj + Ipljj) + (R^ + Ig) 

W2< o4 ) = 04^(R^ + Ig) + 2 04 (R^Rg + I^Ig) + R^ + I^) 

Explicit solution of ^ from equation (49) utilizing the notation defined 
in (50) is: 

^ a(A^XR-x,) 

Equation (51) permits determination of the jg coordinate of parameter plane 
points meeting a bandwidth specification given the value of the oC coordin- 
ate. A negative value under the radical sign infers that no point ^3 

will satisfy the bandwidth specification for the value of assigned. 

w 

When bandwidth is specified at a frequency of w = s/ ^ the bandwidth loci 

in the parameter plane degenerate from conic sections into straight lines 

(for the linear coefficient case). Whether interpretation of these degener 

ate curves is consistent with interpretation of bandwidth curves specified 

at other frequencies has not been established; hence, the reader is caution 

w 

ed to avoid specifying bandwidth at the frequency of 8 / 2 * 

Bandwidth Curves for S-Plane Transfer Functions 

The preceding bandwidth theory is applicable to closed loop transfer 
functions described in terms of the complex variable s. In fact only the 
definitions of the R's and I's need be altered to make z-plane bandwidth 
equations valid for s-plane use. For the sake of completeness the deriva- 
tion below duplicates some of the work in the preceding section. Assme a 
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I 




closed loop transfer function: 



C(S) 



_ K.Q 



y^(o<t^+ (5c,t =t/Sd*4-eJs'' 

CfS) £!(°^Pk*- °<^3« [T«(-C)4- 3/l-g^U|,(~S)] 

■■«« ■ [Tk(-5)+ 

ICsO 

Since bandwidth Is always specified at some real frequency only zeta values 
of zero are relevant. For = 0 the Tj^(O) and Uj^(O) are expressible as: 



T, (0) * cos JlK. 
k' g 



U, (0) = sin 

k g 



This observation allows definitions of the R*s and I’s In equations (47) to 
be altered for s -plane use. 

.K . 



K~0 
K»0 

Rg(Wn> 






2kH 



K-O 

/V 



Ig(» 






£K+/ (52) 



ekfi 



K*o 

N 



IgCWn) -^(~0 ^PK+I 

K*0 



Pk-t-i 



Omitted R's and I's are analogously defined. 

S-plane bandwidth specifications require the magnitude of output to 
Input to have suffered an attenuation A at the bandwidth frequency w, 

D .W* 



Functionally: 



A « 



C(3W.J| 

RUH.)\ 



Evaluate equations (52) at w = w, and substitute the result Into the 
^ n b.w. 

equation above to obtain: 

lid + (R„+ 3X,)| 



A = 






( 53 ) 



From this point all equations (49) and subsequent for z-plane bandwidth 
determination apply Identically for s-plane bandwidth determination. 
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Periodicity of Bandwidth Problem 



A periodic frequency spectrum is characteristic of any signal described 
by z- transforms. Frequently this periodicity characteristic is not intrin- 
sic to the signal being described but is artifically introduced by a fictit- 
ious sampler to facilitate mathematical analysis. The commonplace situa- 
tion for sampled data systems is Illustrated in figure (7) . 




Figure 7 

This frequency spectrum periodicity present in z-transfer functions radi- 
cally alters frequency response interpretation and invalidates or restricts 
prevalent frequency response criteria when applied to pulsed transfer func- 
tions. Among criteria so altered is bandwidth; the 3 db attenuation us- 
ually specified for servo-system bandwidth may never occur when testing 
the pulsed transfer function yet the system itself may be perfectly satis- 
factoiry. This dilemma can profitably be bypassed by return to system des- 
cription in the s-plane even though such description involves transcenden- 
tal functions. Here bandwidth still has meaning and attention need not be 
restricted to sampled signals, a significant advantage when examining the 
system output. 

Sampling makes the previous definition of bandwidth somewhat inade- 
quate (as will subsequently be shown) and agreement on a suitable defini- 
tion of bandwidth is incumbent. Consider the error sampled system in 
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figure (8) . 




Figure 8 

Denote by superscript X*(s) the LaPlace transform of a sampled function 
whose continuous transform is X(s). Reasonable alternatives for definition 
of bandwidth are; 



lc*(s)l . |C(S)| |C(S)| 



The first of these alternatives Is identical to results obtained using z- 
transfer functions and its inadequacy has already been discussed. The 
second of these alternatives allows a different bandwidth for inputs 
and R 2 (s) having the same value at the sampling Instants. For the error 
sampled system shown the output is unaffected by the value of the input 
at other than the sampling instants and one suspects that any definition 
of bandwidth should also have this property. Therefore bandwidth appears 
most suitably defined when stated as: 



. IC(s^ 

Even this statement needs minor modification to be intuitively satisfying. 
The example in the following section will demonstrate a more satisfactory 
and inherently pleasing definition for bandwidth is: 



I lcc$)l 



( 54 ) 
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It is acknowledged that this definition for bandwidth is subject to debate 
when a system is not input or error sampled. 

Block diagram manipulation of the system in figure (8) shows the in- 
put to be related by 



C(s) = 



G(s) 

/4-GW^CS) 



R"(S) 



This relation allows expression of system bandwidth as: 



A 




G(S) 

I+- GH^(S) 



(55) 



A proper approach toward rendering this bandwidth expression fit for para- 
meter plane use might be to examine its individual terms. G(s) is the 
system's LaPlace transfer function which, with the possible exception of 

a hold element, is presumed linear. As such G(s) must be a ratio of two 

sT 

polynomials in "s" and, if a hold element is included, e . GH*(s) is the 
LaPlace transform of a sampled function and it can be represented as an in- 
finite series. 



X*(s) 



X(^r)e 

A^ = 0 



-/vTs 



where x(t) is the impulse response of component X(s). Whenever G(s) and 

H(s) are linear except for hold circuits, as is assumed, the series for 

GH*(s) has a closed form representation which is a rational function of 
sT 

e . Significantly this closed form is readily obtainable from z-trans- 
form tables. 



GH*(s) = GH(z)| sT (56) 

'z = e 

Also noteworthy is that specification of bandwidth at a frequency of s 

jw = jw /2 forces GH*(s) to be real. (i.e. Specify an input-output 
O • w • s 
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attenuation of magnitude A at w /2 in preference to specification of an 

3 

attenuation A' at some other frequency). Such specification reduces the 
problem of separating the real and Imaginary parts of the function 
G(S) 

/-f. &H*(S) 

preparatory to obtaining its magnitude, to a problem of separating the real 
and Imaginary parts of G(jw) which is a rational function of w. 

Still requisite is an equation which allows ready determination of para- 
meter plane points satisfying the bandwidth specification of equation (55) . 
Assume that the linear or product appearance of oC and ^ in the closed loop 
transfer function continues when bandwidth is as specified in equation (55). 
To elaborate, grant that systems having closed loop transfer functions 
shown in equation (46) will retain this ^ame type oC~ jS dependence under 
equation (55) . 

C(S) _ o*:- (-S) -t-/^ (S) 

R^(S) ~ oife(S)4-/3f^,(S)4-oc^fo(S) + f£(5) 

Evaluate the functions f at the specified bandwidth frequency s = jw „ 

D • W • 

and define: 

Rp * Real Part 

Ip = Imaginary Part [fp(jWg y^)] (57) 

Rp- Real Part [fgCjWg^^)] 

Ip » Imaginary Part Ifg(jWp ^ )] 

All R's and I's omitted above are similarly defined. Also define: 

A' - A*T (57a) 

Parameter plane bandwidth loci for bandwidth as defined in (55) can be ob- 
tained using definitions (57) and equations (49) and subsequent which were 
originally developed for z-plane operations. 
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Illustration of Auxiliary Parameter Plane Curves 



The generalization that examples always serve to clarify theory Is 
presumptuous especially in this Instance when bandwidth and steady state 
error loci are already known to be conic sections and straight lines respec- 
tively. Nevertheless an example is presented and surprisingly the general- 
ization holds. In fact a great deal can be learned from this example if 
digression to include several closed loop frequency response curves is 



allowed. Consider the third order sampled data system illustrated in 



figure (2) and direct attention to its closed loop frequency curves 
for representative parameter settings of 



C(314J) 



R*(3W) 



a) 


k = 


0.5, 


kj. - 1.0 


figure (8a) 


b) 


k = 


0.1 


k, » 1.0 


figure (8b) 


These curves 


are drawn 


in figure (8). 


Immediately obvious and disturbing 



is the fact that the closed loop d.c. (w=0) response of this feedback 
system has a magnitude differing from unity. One intuitively expects event- 
ual equality of output and input when a type I system is disturbed by a d.c, 
signal; an expectation which equates to a frequency response of unit magni- 
tude for M-O. The cause for fallacy of intuition and a definition modifica- 
tion which results in more pleasing (familiar) frequency response magnitude 
curves are obtainable from examination of the relative magnitudes of R(s) 
and R*(s). To expand, the magnitude of the frequency response as defined 
for continuous systems and also the definition on which intuition is based 
is 



M = 






R(du») 

whereas a wiser choice for sampled data systems seems to be 



M' = 



C(1UJ) 

R^(JUJ) 
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FREQUENCY RESPONSE CURVES 




FREQUENCY RESPONSE CURVES 
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Now suppose magnitude R(jw) does not equal magnitude R*(jw). Then correct 
intuitive reasoning precludes the same result for sampled data systems. 
However a scale factor may be available to force sampled data frequency 
response curves to rdsemble their counterparts for continuous systems more 
closely. To investigate these conjectures consider the general harmonic in- 
put and form the ratio — 

R»"(S) 



i i KO] = R(s) = 



I 



iirWj R (S) I_q-(S-jujJT 

R(S) ; 

The magnitude of this ratio is of interest at the frequency of the forcing 
function only. 







/ - c 






5 



S = 



= T 



Thus an impulse sampler attenuates a continuous harmonic input signal by a 
factor of »'T”. Such attenuation should be reflected in the frequency re- 
sponse curves; a fact witnessed by the example curves (solid. iTines) of 
figures (8a) and (8b) where T = 0.5. Since attenuation factor "T” is in- 
dependent of forcing frequency w^, dividing the input-output ratio by "T” 
would bring the zero frequency response to unity and the sampled data 
frequency response in general would closely correspond to that for continu 
ous systems. These desirable ends supply the motive for the bandwidth 
definition of equation (55) . Observe that the frequency response curves 
as defined for z-transfer functions 



A = 


C(z) 


_ 


C^(s) 




R(2) 


z = 


/?^Cs) 
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[shovm dashed in figures (8a) and (8b)] have this attenuation factor "T” pre 
sent in both numerator and denominator. Hence these curves are conceptual- 
ly satisfying without recourse to scaling. The curves of figures (9) are 
like those of figures (8) except that s-plane frequency response definitions 
are scaled to have unity value for inputs of = 0 and Y-axis coordinates 
are in decibels. The parameter settings for two of the three curves shown 



a) k = 0.44, 


kj. = 0.177 


figure (9a) 


b) k = 0.39, 


k^ = 1.643 


figure (9b) 



lie on the same parameter plane bandwidth curve in figure (10); consequent- 
ly each frequency curve should exhibit a magnitude of 0.3 (-10.5 db) at 
frequency w = . Although the frequency response curves of figures (9) 

apply to but three specific examples the reader should recognize the extreme 
low pass properties of the hold element and the close correspondence between 



1 


C(S) 




c^cs) 


T 


R^(S) 


and 


R*^CS) 



for w <C which these curves Illustrate. These observations are general- 

ly valid although the frequency at which the generalizations cease to apply 
is dependent on the particular system and the sampling rate. 

These digressions are relevant to parameter plane study to the extent 
that they allow meaningful specification and interpretation of bandwidth 
loci thereon. Several such loci are shown (solid curves) in figure (10). 

For the specified frequency each curve divides the parameter plane into 
regions having attenuation greater than or less than the attenuation evidenced 
by the parameter settings on the curve. If the curve itself represents 
points just meeting a bandwidth requirement then parameter settings corres- 
ponding to points in one of these regions will all exceed the bandwidth 
specification. Moreover no intersection of bandwidth loci occurs for various 
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FREQUENCY RESPONSE CURVES 
FOR A SAMPLED DATA SYSTEM 
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FIGURE 9c 

FREQUENCY RESPONSE CURVES 
FOR A SAMPLED DATA SYSTEM 



Th.e bandwld'th and steady state error loci displayed in figure 10 
are for the third order sampled system shown in figure 2* Equation 
55 governing these bandwidth loci is : 

(1-e-sT) 






3*^ (s+1) (s + 2) 



1 + («056«< ^->31./3 + (.l6o( «.12/3 )e^'^ -W. 026 cK >.1Q/? ) 

e^sT _i,97e2sT^1^20e®T - 0 , 2 '^ 
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magnitude specifications at a given frequency although curves representing 
the same magnitude attenuation specified at various frequencies can conceiv- 
ably have common points. 

The dashed straight line curves in figure (10) represent steady state 
error loci. Perhaps their most significant feature is the vertical slope 
from which the reader infers that steady state error is Independent of the 
parameter /S « 
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CHAPTER V 



Concepts of Dominant Mode Design 

The preponderance of system criteria in the s-plane are based on 
dominant mode analysis and design* Such concepts as damping factor, reso- 
nant frequency, bandwidth and maximum overshoot are meaningful only within 
the context of a dominant mode system. (i.e. What meaning has bandwidth 
if an attenuation of -3 db occurs at three points on the frequency res- 
ponce curve?) Rudiments of dominant mode design reveal the direction any 
extensions of this method must take. 

1. An ideal model must be defined. (For s-plane use, the no- 

zero, two complex pole model is universally accepted as 
standard). ’ 

2. The actual system must be compensated to behave like the 
ideal model. (In the s-plane the actual system approximates 
the ideal second order model if a pair of complex poles 
dominate its behavior). 

Dominant mode design in the parameter plane follows the indicated path. As 
an initial step a z-plane **ideal model" is defined and system behavior (rise 
time, maximum overshoot) is related to model parameters! The following 
section discusses conditions requisite on other systems to insure dominant 
mode (or model^like) behavior. Finally the parameter plane is introduced 
as a vehicle for designing or compensating the actual system to meet the 
requirements for guaranteed dominant mode performance. 

The Dominant Mode Sampled Data System Model and Its Specification 

Choice of a model is critical because dominant mode design can never in- 
tentionally produce systems whose properties exceed those of the model. (To 

illustrate consider an s-plane design in which a fourth order system is 
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compensated to have two dominant complex poles. It may be that the fourth 
order system when compensated outperforms the second order model but such 
performance is coincidental rather than designed). An elaborate model might 
be justified were it not for two additional requirements. 

1. Common specifications must be readily translatable as model 
parameters (root locations). By way of example, a maximum 
overshoot specification can be translated into a zeta re- 
striction on the roots of the second order s-plane model. 

2. The model's pole-zero configuration must have simplicity 
sufficient to enable the designer to recognize when model 
like behavior (dominance) is achieved. Suppose, for an 
example by contrast, it is possible to determine a fourth 
order ideal model system. A designer now compensates an 
existing system to approach this model; but approach it 
where? Should he design to have two poles in close prox- 
imity ahd. two further removed or should all poles be moder- 
ately removed from those of the model or is some other ap- 
proximation to be preferred? 

These considerations indicate an acceptable model to be one whose closed 
loop z-transfer function is 



C(g) ^ k(z--z,) 

R(z) (Z-ff)(z-P,) 



where z Is a real zero 
1 

fil'd are conjugate complex poles 

- // ~P i )( l-Pi) This choice of "k" guarantees a type I system 
\ / “"Z|/ 

with zero steady state error to unit step inputs. 
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The utility of this model is contingent upon the facility with which 
system specifications are translatable into acceptable pole-zero configura- 
tions. To amplify, unless acceptable model pole-zero configurations can 
be simply determined design mights prof itably proceed in some direction other 
than in forcing the system to resemble a 'model. The following paragraphs 
develop the required relations between maximum overshoot, rise time and 
model pole-zero locations. This development largely duplicates the work 
on pages 170 through 179 in Kuo and is included to accentuate approximations 
inherent in the model itself [2]. These approximations are over and above 
those made when compensating higher order systems to a dominant mode and are 
not required for continuous systems. 

The rise time and maximum overshoot of the system described by equation 
(58) to a unit step input are desired. This requires expression of c*(t). 

/-7'S — ^ ^ ~Z- 

(7-f^)(Z-p.) ' -i-t 



By the inversion integral: 



C(nT) = 



I f 

STTiJ (Z- l)(2- P,)(z-f,) 
r 



2 






d 



Z 



where P is a circle of radius e centered at the origin and 
poles of C(z). Applying the residue theorem to this integral 



enclosing all 
one obtains 



CCnT) = . K(F. -z.) 

0- ?)(/-!;) (i?-0(f;-^) (f-O(p-R) 

which under definition of 



p, = lP| 

becomes 

c(nT) = ;_ 
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or 



cos -©i + n § 



c(nT) = I -h 2 k 



\(p-z.) P" I 
I(P-0(P,-7^)| 



where Y = arg(pj^-Zj^) - arg(pj^-l) 



-0j - arg(p^-z^) - arg(p^-l) - ^ 

With specification of the closed loop poles and zero the above expression 
is sufficient to determine the output c*(t) at the sampling instants. How- 
ever this expression has two serious defects. 

1. output is specified only at the sampling instants whereas 
peak overshoot in the continuous output is likely to occur 
during the intersampling period* 

2. The output c*(t) is specified as a function of three para- 
meters, a) real part of complex poles, b) imaginary part of 
complex poles and c) location of real! zero; consequently its 
universal solution is difficult to display. That is to say - 
system criteria which are functions of three variables are 
not conveniently represented in two dimensions or less. 

The next step is circumvention of these difficulties. The first problem 

can be attacked in several ways among which are employment of modified z- 

transforms and approximation of the continuous output c(t) by a continuous 

function c (t) which agrees with c(t) at the sampling instants. The complex- 
a 

ity of the modified z-transform and the fact that zeros of C(z,m) are func- 
tions of m preclude the use of this avenue. Approximation of the continu- 
ous output with a continuous function which agrees at the sampling instants 
offers many advantages. Among them, the peak overshoot of this approximating 
function can be determined and if, as is hoped, the approximation is good 
then this peak should be close to the true peak overshoot i The second pro'y^ 
blem is somewhat more difficulty. It Involves specifying the poles and zeros 
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in such a manner that overshoot can be expressed in terns of two parameters 
i 

instead of three. Offhand there is no direct procedure but rather it is 
one of trial and error. The following examples should clarify what is being 
attempted. 

1. For continuous second order systems the complex poles can be 

expressed as (T+ jw or w^( ^ j ). When the poles 

are expressed by ^ and w^ then overshoot becomes a function 
of ^ alone. 

2 2 

2. Suppose ~ oi.^+ Y • It is obvious 

that £ can be expressed as a function of two variables by de- 
fining ^ . Then 

- S^-S + Y 

For the second order sampled data system of interest a successful choice of 
variables has been found. Define: 

- 2 ,) - ^ 

Expression of c(nT) in terms of these parameters is most easily accomplished 
by recourse to geometry. 

Im[z] 



Re[z] 

Figure 11 

Consider figure 11. By Inspection the folldwing definitions can be made. 
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Hi = argument (p^-1) 
argument 

$ = argument (p^) 

arg - arg (p^-1) + y = % ^ 

These definitions of J and ^ are consistent with those made earlier. 

An Interesting representation for the secant of ^ results from these 
definitions. Define the secant in terms of the legs of right triangle (p^, 
z^, b). 



sec ^ = 






Z,-b 



Now observe that the area of triangle (p^, z^, 1) can be expressed in two 
ways 



area A (p^, z^,l) 




Zl-b|=^ \P,-P,\ \ I-Z, \ 

Therefore 




-i(|z,-bl|i-R|) 



sec 









- e 1' 




iP-z,| 




1 |-Z,| 'l'^- 




= S k 











But this expression is identical to that which appears in the equation for 
c(nT). Hence c(nT) can be written as: 

c(nT)=|seC^| |p, l^cos (a' $+ f'-ir) 

The previously defined representation of "Pj^" as a function of ^ and w^ 
allows expression of c(nT) as: 

c(nT) =|sec ^ I cos(i//-^^ (4, nT 

At this point it becomes a simple matter to fit a continuous curve through 
the discrete output c(nT). To do this simply substitute the continuous 
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variable t for the discrete variable nT. 



,(t) = Isec ^1 cosif^ bJA ^- t) 



( 59 ) 



c^(t) Is the continuous function which will be used in place of the true 

output c(t). While c (t) coincides with c(t) at sampling instants there 

a 

is no requirement that it coincide during the intersampling periods as 

well. Nevertheless it is assumed that c^(t) is "close to" c(t) at all 

times and that the system output is adequately represented by c^(t). Under 

this assumption c^(t) can be differentiated to obtain peak time, maximum 

overshoot and other characteristics inherent in c(t). 

The peak time is found by differentiating c (t) with respect to t and 

a 



equating the resulting expression to zero. 



I 



m2ix 



ui. 









(60) 



The maximum overshoot is found by substituting the expression for T into 

• max 



c (t). 
a 



(»ax) Isec si e ^ '■’il/ 

I u — — — 



-h / 



(61) 



The portion of this expression underlined is the maximum overshoot to a 
unit step input. Charts of overshoot and — • as functions of ^ and 
C can be found on pages 175 and 176 of Kuo [2], 

The equations for T (60) and c (max), (61) characterize a second 

max a 

order sampled data system in terms of two parameters. However they are not 
exact representations and the problem is not really complete until regions 
of validity for this characterization are determined. Lindorff (by what 
seems to be Intuitive argument) has determined that the approximation is 
good if the complex poles are in the first and fourth quadrants and are not 
too close to the origin of the z-plane [9]. [Poles in the first and fourth 
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quadrants Insure at least four samples per cycle of the dominant oscilla- 
tory mode and thus give a reasonably good approximation for this mode. How- 
ever, if the poles are too close to the origin the rapid decay offered by 
the |p|” factor will predominate to the point where the z-transform [and thus 

the approximating function c (t)] no longer adequately characterize the 

a 

system. That is, the sampling rate is too slow to follow the decay]. 



Analysis of Dominance Requirements 

For the z-plane dominant mode model developed in Kuo and reviewed in the 
preceding section to have value the cpnditions under which model-like be- 
havior prevails must be determined. These requirements can be established 
by considering the Inverse z-transform integral. 



c(nT) = 



rf?/ 



CCz)-z"-' dz 



where P is a path which encloses all singularities of C(z). But 



C(z) * G(z)R(z) 

C(3) - ^ 

tn-F,) z-i 

Substituting this expression into the integral equation yields 



for a unit step input. 



c(nT) 



I 






II (Z-Zk) 

Ji=! 



z-l 



:“d: 



II (z-fl) 

This integral can now be evaluated by the theory of residues. Expressed in 



this manner it becomes: 



II (f^ ~Zk) 



1 state output 

IT(^-py) ilO-O' 



Analysis of this expression reveals the two conditions necessary for dominant 
mode behavior. 

1. All but the two dominant poles must lie near the origin of 
the z-plane. After a small number of samples the residues 
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contfibuted ;>y poles near the origin are negligible; the p” 
term dominates these residues and as a factor it becomes 
almost zero. 

2. Excess zeros over the real zero employed in the second order 
model must lie near the origin of the z-plane. When so 
located two effects are assumed to result. 

a. The residues of poles near the origin (which are 
neglected in the second order model) will be small. 

b. These zeros form dipoles (pole-zero pairs close to- 

gether) which have almost negligible effect on the 
residue of the dominant poles, [i.e. If ~ 

(p^-Pj^) then these common factors cancel in the ex- 
pression for the dominant pole (Where time de- 

lays are absent, dominant mode design requires as 
many poles as zeros to be ignored or lie near the 
origin. This proximity of ignored poles and ignored 
zeros essentially" guarantees dipole formation when 
viewed from the location of the dominant poles) . 

Principles of Dominant Mode Compensation in the Parameter Plane 

With model behavior and restrictions to insure model behavior in higher 
order system established attention turns to the problem of properly compensat- 
ing existing systems to resemble the model. Parameter plane methods are well 
suited to this aspect of dominant modeidesign. This approach not only allows 
selection of the dominant (model) poles but also guarantees dominant mode (or 
model like) behavior by forcing compliance with the restrictions mentioned 
above. Tacit is the assumption that parameter values ( o(. , /^ ) exist which 
allow accomplishment of this ambitious objective. Justification of these 
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statements can be had by recourse to Chapter III. Recall that the image 
of a 2-plane contour divides the parameter plane into regions whose distin- 
guishing characteristic is that any parameter setting within a given region 
places an equal number of characteristic equation roots within the z-plane 
contour which was mapped. Suppose a circular contour centered at the origin 
of the z-plane and having a radius small enough so that all poles within 
it can be considered negligible, is mapped through the characteristic equa- 
tion into the parameter plane. Of the parameter plane regions defined by 
this mapping the dominant mode designer considers only those for which all 
roots but two of the characteristic equation lie within the z-plane contour. 
If within such a region parameter settings exist which place two poles near 
the unit circle then these poles must be dominant (because it is known that 
no other poles except those near the origin exist) . To find these poles 
auxiliary contours of constant w and constant C are mapped through the 
characteristic equation Into the parameter plane. The Intersection of an 
acceptable contour and $ contour in the region mentioned above constitu- 
tes a set of parameter values which simultaneously guarantees dominance and 
the desired pole locations. 

To this point the closed loop zero specification required by the domin- 
ant mode model has been' ignored. Such slight is intentional since zeros of 
the open and closed loop transfer function are identical and are not general- 
ly subject to modification by available system parameters. Therefore dis- 
counting cancellation, zeros should be considered as entities which allow or 
disallow dominant mode compensation rather than as entities to be modified. 

Dominant mode design using parameter plane concepts is readily extended 
to continuous systems described in the s-plane. Instead of a small circle 
near the origin a contour of constant (T ( S“ =Cv ) sufficiently removed from 
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the imaginary axis is mapped through the characteristic equation. The image 
of this s-plane contour divides the parameter plane in exactly the same 
manner as the image of the small circle in the z-plane; subsequent s and 
z-plane procedures are identical. In fact dominant mode design in the s- 
plane is generally simpler than in the z-plane because the s-plane dominant 
mode model has no zeros in the numerator of its transfer function. 

Cascade Compensation (for Dominant Mode Design ) 

The principles discussed in the preceding chapter fairly well dictate 
the form that a cascade compensator must have if dominant mode design is 
to be achieved. Consider the sampled data system of figure (12). 




G(z) is the plant transfer function which is assumed to be of the form; 

(62) 

IT (z + 0 

The system open loop zeros provide the best clue as to the feasibility of ' 
dominant mode design. Since the zeros of the open and closed loop transfer 
functions are identical, all zeros except one must lie near the origin of 
the z-plane. When this condition is not met no cascade compensation is 
capable of driving the system to the second order dominant mode desired. 

It is also prudent to observe the system poles before undertaking dominant 
mode design. Cascade compensation can succeed, disallowing pole-zero can- 
cellation, only if addition of poles and zeros offers the possibility of 



G(z) = 



80 



reshaping the root locus so that all but two complex poles and one real zero 
lie near the origin. If the general uncompensated system pole-zero configura- 
tion is such that this is manifestly impossible then there is no point in 
employing parameter plane procedures to prove it. 

These observations can be made more rigorous by translating them into 
restriction on^the transfer function of the cascade compensator. 

1. The zeros of the cascade compensator must lie near the origin 
of the z-plane (assuming that one system zero is not so 
located). In this way these zeros can be neglected without 
affecting the validity of the dominant second order approxi- 
mation. 

2. The numerator should be of the same order in z as the de- 
nominator (primarily because time delay is an undesirable 
property in a compensator) 

3. The steady state gain of the compensator should be unity. 

With these properties in mind the transfer function of the cascade compen- 
sator must have the form: 



_ • 4“ Q, -h ^ 



7T+ Z 



(63) 



( 2 | ~Z. G.Q 

With this cascade compensator installed the closed loop transfer function 
of the system in figure (12) is expressible as: 



C(z) ^ 
R(2) 



M 

k 4- - - - -I- a, 4- a®) 



+Qo) / I (z-pj + ■ +ao)2'" /I (z-z,) 



(64) 



x=i 



J- 



where k, z . , and p are known fixed parameters of’jfthe plant. The a 's are 
i i » j 

now chosen to yield a satisfactory dominant mode second order system. While 
equations (63) and (64) appear formidable, parameter plane methods will yield 
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a solution almost by inspection especially if but two of the 3j*s are 
variable. 

This method is also applicable when pole-zero cancellation is allowed. 
In this case the cascade compensator has a transfer function of the form: 

Q (z) = -ha, -h ao)z'^ 

x=i 

where (z + r^) and (z + t^) cancel poles and zeros of the plant transfer 
function respectively. 



Feedback Compensation (for Dominant Mode Design ) 

Feedback compensation instead of cascade compensation can be used to 
effect dominant mode design. As with cascade compensation certain restric- 
tions are placed on the feedback transfer function if the chosen second 
order dominant model is to be valid. Consider the sampled data system of 
figure (13) . 




where h” » numerator of H(z) 

= denominator of H(z) 
g” * numerator of G(z) 

G^ - denominator of G(Z) 



Figure (13) 

Assume that: 

a. The plant transfer function G(s) contains at least one integra- 

tion. (l.e. The denominator of G(s) contains a factor of s 
or the denominator of G(z) contains a factor [z-l[). 

b. At most all but one zero of the plant transfer function is 

located near the origin of the z-plane. 
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Requirements: 

1. Design a ^feedback compensator H(z) which yields a second order 
dominant mode system with acceptable time response characteristics. 

2. The system must follow a step input with zero steady state error. 
The first step toward solution will be to apply the assumptions with the 
objective of obtaining an acceptable class of transfer functions H(z). 
Parameter plane analysis applied to this class will then yield acceptable 
H(z)'s if they exist. The closed loop transfer function for the system of 
figure (13) is 



C(Z) 6(Z) 

R(Z) / + G(z)HCz) 

C(Z) 

R(z) 

The first restriction on H(z) can be obtained by observing that poles of 
H(z) are zeros of the closed loop transfer function. To effect dominant 
mode design these closed loop zeros [or equivalently poles of H(z)] must lie 
near the origin. The second restriction on H(z) results from the requirement 
that there be zero steady state error to a step input. This restriction is 
obtained by applying the final value theorem to the closed loop transfer 
function. 



c( oO ) 

c( oO ) 



r :z- 1 ) 1 ^( 2 ) n = 



/ for unit step 

input 



lif^ 

Z-^l 



r (2-1) z ~ i ^ r &'*><■' ~ ) 

[_ 2 (z-O GW+G^J iZi^JFTGVPi] 



1. iiD d 

But z-*i [G ] ^ 0 since the plant transfer function was assumed to contain an 
integration. Hence: 



c( oO) 



IV^' 

Z-*l 



j/M 

z-»; 
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or 



lim 



H(z) = 1 



To meet these requirements the feedback path must have a transfer function 
of the form m 

(Z-')E 

H(z) - 1 + (65) 

This transfer function Is slightly more restrictive' than necessary since all 
poles of H(z) were placed at the origin (for computational ease) rather 
than near the origin which Is the specific requirement. Also notice that 
H(z) is written as containing no time delays although this is not a specific 
requirement. Time delays can easily be added by Increasing the powers of 
z in the denominator of H(z) . 



Example of Dominant Mode Design in the Parameter Plane 

t 

A numerical example is appropriate to Illustrate theory and to rein- 
force stated results. Suppose the plant shown in figure 14 is to be compen- 
sated to a dominant mode which satisfies the specifications 

1.2 <" M ^ C 1.4 

rise time - less than 5 sampling periods 
As an Initial step in the compensation procedures the zeros of th$s plant's 
closed loop transfer function (equation 66) are checked for compliance with 
necessary conditions for dominant mode design. Namely, validity of dominant 
mode design presupposes all zeros but one are sufficiently near the origin 
to be negligible. The real zero at z » -.295 stretches this requirement but 
when compared to the model zero at z = -3.72 it will probably not be suf- 
ficient to destroy dominant mode behavior. Tables describing model be- 
havior such as those referenced in the first section of this chapter are now 
consulted for ^ values which satisfy the specifications. For this 
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/ 

FIGURE 14. ^ 



BLOCK DIAGRAM OP A FOURTH ORDER 
PLANT AND PROSPECTIVE COMPENSATION 
EMPLOYED TO EFFECT DESIRED DOMINANT 
MODE SYSTEM BEHAVIOR 



The closed loop transfer function of the system shown above 
is: 

C(z) _ N(z) 

" D(z) 

where 

N(z) = 0 . 022 (z 2 )(z-i- 3 * 72 )(z+ 0 . 29 ^)(z+ 0 . 015 ) 

A c' (6 

D(z) ss z^-l- (0.0220/8 *1.814.) z 5 

(0. 0220o^ +0.0667/?+!. 066) z^ 

+ (0.o 667<< -0.0632/* -0.116) z 3 

+ (-0.0632 o( -0.0252/3 4-0, 00128 )z2 

f 

+ (-0.0252 -0.00037/3 )z 

+ (-0.00037 ) 



\ 

\ 
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example approximate value ranges 
-30* < S < 10“ 

0.3 ^ ^ < 0.45 

appear satisfactory. 

The parameter plane is now employed to effect dominant mode behavior in 
the actual system. Those system poles whose settling time is five times or 
more faster than that of the dominant mode will be considered negligible. 
Since w (the settling time Indicator) of the dominant mode is about unity 
all other poles must lie within a circle of radius w =>0.2 centered at 
the origin. The image of this settling time contour divides the parameter 
plane into regions containing an equal ninnber of roots to the characteris- 
tic equation (figure 15); the region for which all but two characteristic 
roots lie within the designated settling time contour is shaded. The pre- 
sence of constant ^ curves in this region within the range needed to meet 
specifications is ascertained and on each curve the range of satisfactory ^ 
values encountered in traversing the region is marked. These actions serve 
to delineate a parameter plane region (shaded in figure 16) within which all 
parameter settings provide a second order dominant system meeting specifica- 
tions. Figure 17 compares the transient response of a true second order 
system with the response of the system in figure 14 compensated for dominant 
mode behavior. The dominant poles and zeros in the latter case identify 
with the poles and zero of the true second order system whose response is 
displayed. Numerical values for the compensated system are as follows: 
parameter settings 

oC » -0.13 /S » 0.73 

dominant closed loop poles and zero 

poles: z = 0.776 + jO.362 zero: = z = -3.72 
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ignored closed loop poles and zeros 
poles*- z ■ 0.182 + j0.078 
z » -0.016 
z » -0.106 



zeros: z * -0.295 
z » -0.015 
z^ “ 0.0 
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CHAPTER VI 



Frequently the niceties of a compensation procedure are more clearly revealed 
by numerical examples. In deference to this hypothesis several examples 
emphasizing various aspects of parameter plane theory are presented in this 
chapter. 

Example of Compensation by Parameter Plane Technique 

The fragmentary examples in preceding chapters relating to the system 
shown in figure (2) can be gathered into a fairly detailed and representative 
illustration of parameter plane technique. Typically the engineer is given 
a basic plant such as that shown in figure (2) and he is told-, to modify its 

it . . 

performance to meet specification. Suppose these specifications '^require 
maximum overshoot to step inputs to be under twenty per cent, limit the 
sampling rate to two samples per second, and require steady state errors 
to unit ramp inputs of less than one. In addition short rise time and short 
settling time is desirable. Direct consideration of all these specifications 
is not practical in the parameter plane but tractable indicators for each 
exist. These specifications modified for parameter plane utilization be- 
come: 

> 0-5 lieu of < 1.2 

2. E cT 1 direct consideration 

s.s. 

3. large bandwidth in lieu of short rise time 

4. short settllnt time direct consideration 

The engineer decides that if feasible a combination of gain adjustment and 
derivative feedback offers the most attractive way to meet these specifica- 
tions. The closed loop z-transfer function with the perspective compensa- 
tion Incorporated (equation 32) is amenable to parameter plane analysis. 
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Pertinent curves .Sufficient to determine the parameter settings required to 
effect a satisfactory system are shown in figure 18. One observes that to 
meet overshoot specifications parameter settings ( oC , ) must lie to the 

left of the 5 “ 0.5 contour while steady state error requires these set- 
tings to be to the right of the error locus shown. Unfortunately the inter- 
section of these two areas is rather small leaving the designer little room 
to maneuver for the desirable properties of short settling time and short 
rise time. Nevertheless a constant settling time curve is shown; all para- 
meter settings within the contour yield a system with a shorter settling 
time than is obtained for parameter settings outside the contour. Band- 
width curves which are Indicative of rise time are also drawn; regions to 
the right of each curve have shorter rise times than those to the left. 

From these parameter plane considerations one concludes that settings of 
o6= k = 0.5 
/S “ kk^ "0.5 

should yield an acceptable system. The time response (at the sampling 
Instants) for this choice of parameters is displayed for step and ramp in- 
puts in figure 19a and 19b respectively. 

Example of Analysis by Parameter Plane Technique 

In addition to its capabilities as a design tool the parameter plane 
also .provides an excellent vehicle for analysis. Consider the second and 
third order feedback systems shown in figure 20a and 20 b respectively and 
suppose them to be components of a still larger system. If pole b is far 
removed from pole a along the negative real axis it is well known that the 
third order system will behave essentially like its second order counterpart. 
However approximate second order behavior of this component depends not only 
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on the distance that pole b is "removed" but additionally it depends on 
gain values k and k^. A knowledge of when second order behavior prevails 
as a function of k, k^., and b is of consequence in analyzing the overall 
system. Furthermore, if k and k^ are used to adjust overall system perfor- 
mance - for what range of values can these adjustments be treated as ad- 
justments to a second order system and when must the third pole be consider- 
ed? These questions are conveniently answered on the parameter plane, 
shown in figure 21, where it is assumed that nominally these poles ate ad- 
justed to a zeta value of 0.5. (If additional ^ values are of interest 
these contours are also mapped into the parameter plane) . Proximity of 
second and third order C = 0-5 curves denotes approximately equal charac- 
teristic equation roots for parameter settings in these regions; consequent 
ly the third order system is closely approximated by its second order count 
erpart. Figure 21 shows such approximations to be valid for 



0.6 < oC 4 k <1.6 
0 < jS ^ kk^< 0.6 



for b = 10 



0 < f3 <0.35 

Notice that aven a 10:1 ratio between poles a and b does not justify second 
order approximation torci> 2.4 or ^ > 1.5. In figure 22 the effect of 

these approximations on bandwidth is shown. 

The author cautions against haphazard extension of these approxima- 
tions. For instance, proximity of ^ - curves for a second and fourth 
order system would not be sufficient to guarantee the validity of a second 
order approximation. Not only must the location of the other two roots to 
the fourth order system be investigated but the region of proximity for the 
^ - curves must be checked for a corresponding proximity of w^ values as 



well. 
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Advanced Example of Compensation by Parameter Plane Technique 

As mentioned In the introduction most graphical methods are seriously 
limited when called upon to display system performance measures as a func- 
tion of more than one variable. In this respect parameter plane methods 
have definite advantages for performance indicators such as damping, settling 
time and bandwidth are Inherently displayed as functions of two system vari- 
ables. Furthermore, when only a single zeta value (or a single settling 
time or a single bandwidth) is of interest parameter plane methods compet- 
ently display the effect of three variables on system performance. Even 
more variables can be considered without undue effort as will be illustrat- 
ed in designing compensation for the system of figure 23. Here the design- 
er has four parameters at his disposal, the system gain, the amount of de- 
rivative feedback, and the pole and zero locations of the single section 
filter. Settings for these parameters which allow the system to meet the 

V 

following performance specifications are to be determined. 

1. Maximum overshoot must be less than 1.25 or zeta is to be 
equal to or greater than 0.4. 

2. The steady state error to unit ramp Inputs must be less 
than 0.25. 

3. The bandwidth (-3 db point) of this system must occur at a 
radian frequency greater than 7T to allow sufficiently fast 
response. 

Examination of the closed loop transfer function reveals its suitability 
for parameter plane analysis with oC and yd capable of definition in any 
of several ways. The curves Illustrated result from the choice 
o4 « k, yd =kp/q 

with p and k^ considered as third and fourth parameters. Initially k^. is 
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equated to zero and the efficacy of lag and lead compensation is invet'cl- 
gated. The effect of such compensation on stability and bandwidth, as 
displayed in a parameter plane format, is shown in figures 24 and 25 res- 
pectively. That neither lag nor lead compensation is capable of rendering 
a satisfactory system is immediately evident. The intersection of the para- 
meter plane region denoting stability with that region providing a satis- 
factory steady state error is empty for lead compensation whereas the inter- 
section of those regions denoting stability with those denoting acceptable 
bandwidth is empty for lag compensation. However the parameter plane dis- 
play does indicate that a "good location" for a lead section pole is around 
s “ -8 and that a lag pole might profitably be located at s = -0.05. 

The effect of derivative feedback in combination with lead and lag 
filters can now be investigated. The damping curves of figure 26 and the 
bandwidth curves of figure 27 display the results for various values of de- 
rivative feedback used in conjunction with a lead filter (or more accurately 
a single section filter whose pole is located at s - -8). Figure 28 dis- 
plays an analogous result for feedback used in conjunction with a single 
section filter whose pole is located at s = -0.05. 

In both instances system performance is enhanced by employment of some 
derivative feedback. In this case lag filters provide the easiest way to 
improve steady state accuracy but even with feedback their bandwidth limita- 
tions (slow rise times) cannot be overcome. On the other hand lead filters 
readily provide acceptable bandwidth and with derivative feedback both damp- 
ing and error specifications can be met. Intersection of acceptable damp- 
ing, error, and bandwidth regions indicate parameter values of 



k = oL = 5.0 



kj. = 0.5 



P 



= 8.0 




should suffice 
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CHAPTER VII 



The general availability of the digital computer has fostered revolu- 
tionary advances in servo-system and especially, in sampled data system 
design. However many of these new approaches completely discard rather 
than refurbish existing methods and thus void much of the experience gained 
by practising engineers. In contrast parameter plane methods retain the 
principle of design by such familiar specifications as damping and settl- 
ing time but where classical methods could display the effect on only one 
variable on these specifications, parameter plane methods can display the 
simultaneous effect of two. This paper has extended parameter plane methods 
in two directions. First the requirement that parameters appear linearly 
in the characteristic equation was eased to permit consideration of those 
cases in which their product is also present. Extension in the second 
direction increased the Information which could be displayed, and hence de- 
signed for, on the parameter plane to include bandwidth and steady state 
error specifications. This latter extension is particularly significant 
because it relates performance to the zeros of the closed loop transfer 
function which were previously ignored. Another novel capability of the 
parameter plane which was introduced in Chapter 5 is the simultaneous de- 
sign of a dominant mode system and the guarantee of such dominance. 

A great deal of further investigation into parameter plane methods in 
both theory and application still remains. To amplify questions typical of 
those to be resolved consider the consequences of allowing product terms to 
appear in the coefficients. Now points in the s- or z-plane no longer have 
unique images in the parameter plane and some points in the transform planes 
could conceivably have no real image at all. It is easy enough to postulate 
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that available parameters cannot place characteristic equation roots at 
these latter locations and that closed contours in the transform planes will 
result in consistent contours in the parameter plane. That is - these para- 
meter plane contours are expected to be consistent with the interpretation 
expounded in Chapter 3. While such postulation seems reasonable enough 
for the product case, it certainly warrants a more rigorous framework if 
these methods are to be extended still further. 

To the author's knowledge the bandwidth and error curves represent the 
first display of information on the parameter plane which is not directly 
obtainable from the characteristic equation. Research to discover other 
performance Indicators capable of being displayed on the parameter plane 
should certainly be fruitful. In fact the author has obtained the equa- 
tion governing loci of constant root sensitivity, but each point thereon 
must be computed by an iterative procedure and it is too unwieldy to be 
generally useful. Utilization of this auxiliary information in conjunction 
with that normally available on the parameter plane is Illustrated in 
Chapter 6. Another facet of parameter plane methods illustrated by these 
problems is the relative ease with which the effect of a third variable 
parameter can be analysed. Thiii ease represents a significant improvement 
over the almost unworkable situation a three variable problem presents to 
root locus methods. 

The author feels that parameter plane methods as enhanced by this paper 
should equal or exceed, for two parameter problems, those capabilities which 
the root locus brings to problems of one parameter. Furthermore, the re- 
liance on familiar performance indicators minimizes prerequisite studies 
and enhances the experience of most engineers. These trait^ auger its more 
widespread appearance. 
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APPENDIX 



Some Self-Explanatory Parameter Plane Curves 

These self-explanatory curves should serve to relate the familiar 
second order system to some parameter plane representations of it. 

A knowledge of the general pattern of such relations will serve the 
reader when interpreting parameter plane cuirves for more complex 
systems. 
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